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Abstract 

This  study  examines  reduced-order  Kalman  filters 
designed  to  improve  performance  of  the  F-4E/G  long  range 
air-to-air  missile  capability  (LRI  function).  Operational 
requirements  dictate  a  high  degree  of  accuracy  and  con¬ 
straints  imposed  by  existing  hardware  mandate  minimal  com¬ 
plexity.  Two  linear  dynamics  models  are  proposed,  one  based 
on  constant  target  velocity,  and  the  other  based  on  time- 
correlated  target  acceleration.  Both  are  defined  in  iner¬ 
tial  Cartesian  coordinates  aligned  with  north,  east,  and 
down.  A  nonlinear  model  is  developed  for  measurements 
available  in  the  existing  F-4E/G  hardware,  including  range, 
range  rate,  radar  antenna  gimbal  angles,  and  radar  antenna 
rates.  The  models  are  implemented  in  extended  Kalman  filter 
formulations  employing  linear  propagation  equations  to  avoid 
on-line  numerical  integration.  Performance  evaluations  are 
performed  on  three  test  trajectories  using  Monte  Carlo  anal¬ 
ysis.  Filter  tuning,  error  budgets,  adaptive  techniques, 
and  observability  issues  are  addressed  during  filter  evalua¬ 
tion.  Results  of  the  evaluation  indicate  the  filter  designs 
can  meet  the  requirements  of  the  F-4E/G  fire  control  system. 
Recommendations  are  made  for  continued  testing  and  for  oper¬ 
ational  implementation. 


KALMAN  FILTER  DESIGN  FOR  THE  LONG  RANGE  INTERCEPT 


FUNCTION  OF  THE  F-4E/G  FIRE  CONTROL  SYSTEM 


I.  Introduction 


Background 

Fruitful  applications  of  Kalman  filtering  in  the  area 
of  airborne  target  tracking  are  well  documented  in  current 
technical  literature  (1,4).  Extensive  research  has  been 
conducted  in  this  area  since  the  advent  of  Kalman  filter 
theory  in  the  early  1960s,  and  today,  fully  developed  sys¬ 
tematic  design  procedures  are  available  (1:341).  Neverthe¬ 
less,  engineers  involved  in  specific  implementations  face 
recurring  difficulties  associated  with  the  limited  speed, 
memory,  and  wordlength  of  the  small  on-line  computers  used 
to  implement  real  time  filters.  This  thesis  is  an  attempt 
to  build  on  the  results  of  previous  research  by  applying 
established  Kalman  filtering  techniques  to  a  highly  specific 
problem,  estimating  the  velocity  of  long  range  targets  under 
the  contraints  imposed  by  hardware  in  the  fire  control 
system  of  the  F-4E/G  aircraft. 

The  F-4E/G's  fire  control  system  uses  an  onboard  radar, 
the  APQ-120,  to  track  airborne  targets  prior  to  launching 
air-to-air  missiles  (9,11,21).  A  multifunction  microproces¬ 
sor  controls  the  fire  control  system  through  associated 


software,  designated  the  Operational  Flight  Program  (OFP). 


One  function  of  the  OFP  is  to  combine  radar  target  data  with 


F-4E/G  ownship  attitude  and  velocity  data  to  obtain  accurate 


estimates  of  target  velocity  and  acceleration.  These  esti¬ 


mates  are  then  used  to  determine  when  the  target  is  within 


an  envelope  for  successfully  launching  a  missile.  Cur¬ 


rently,  the  OFP  uses  a  Wiener-Hopf  filter  to  estimate  target 


velocity  and  acceleration  for  use  in  computing  the  maximum 


and  minimum  lethal  ranges  of  a  missile  under  current  flight 


conditions.  The  output  of  the  OFP  then  drives  a  steering 


dot  and  a  maximum  missile  range  indicator  in  the  cockpit. 


These  instruments  provide  the  pilot  with  the  information 


needed  to  set  up  a  successful  attack  against  distant  targets 


(the  long  range  intercept,  or  LRI,  function)  and  closer  air 


combat  maneuvering  targets  as  well  (the  ACM  function). 


Ideally,  the  pilot  should  be  able  to  maneuver  the  aircraft 


by  following  the  steering  dot  and  to  launch  a  missile  suc¬ 


cessfully  when  the  range  indicator  shows  the  target  is 


within  range.  However,  pilots  have  complained  that  the 


steering  dot  becomes  unstable  and  the  maximum  missile  range 


indicator  temporarily  goes  to  zero  whenever  the  F-4E/G  per¬ 


forms  high  roll  rate  combat  maneuvers  (5,8,9). 


The  steering  dot  instability  stems  from  two  sources,  an 


increased  noise  level  in  the  radar's  rate  gyro  outputs  during 


aircraft  maneuvers  and  the  tendency  of  the  Wiener-Hopf  filter 


to  amplify  this  noise  (5:1).  Improvements  might  be  achieved 


by  replacing  either  the  radar  or  the  filter.  However, 
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replacing  the  APQ-120  radar  is  a  costly  and  uncertain  fix  at 
best,  and  no  funds  are  available  for  this  option.  The  second 
alternative,  replacing  the  Wiener-Hopf  filter,  could  be 
accomplished  at  much  lower  cost  by  reprogramming  the  OFF. 

For  this  reason,  it  is  proposed  to  design  a  Kalman  filter  to 
replace  the  current  Wiener-Hopf  filter.  If  such  a  Kalman 
filter  could  be  feasibly  implemented  in  the  constrained 
hardware  of  the  F-4E/G  fire  control  system,  it  could  opti¬ 
mally  process  the  poor  quality  radar  data  in  such  a  way  as 
to  eliminate  the  steering  dot  instability  without  sacri¬ 
ficing  performance. 

Problem  Statement 

The  purpose  of  this  study  is  to  determine  the  feasi¬ 
bility  of  implementing  a  suitable  Kalman  filter  in  the  F-4E/G 
fire  control  system  under  the  constraints  imposed  by  current 
hardware.  It  is  well  established  that  higher  order  Kalman 
algorithms  can  achieve  improved  estimation  performance  over 
filters  of  the  Wiener-Hopf  type,  providing  that  sufficient 
memory  and  computational  speed  are  available  from  the  com¬ 
puter  (4,  6:157).  In  this  respect,  the  F-4E/G's  computer 
imposes  severe  constraints.  It  operates  at  300,000  opera¬ 
tions  per  second  (300  KOPS),  which  is  slow  compared  to 
state-of-the-art  microprocessors  for  real  time  avionics. 

For  instance,  the  fire  control  system  in  the  F-16  currently 
operates  at  700  KOPS  (7).  A  maximum  filter  cycle  time  is 
specified  as  40  msec  by  OFP  software  engineers  (8).  This  is 


actually  a  current  requirement  of  the  Wiener-Hopf  filter 
needed  to  produce  estimates  of  the  desired  accuracy  and 
still  leave  time  for  other  steering  and  range  computations 
required  to  update  the  steering  dot  and  range  indicators  in 
the  cockpit  at  a  10  Hz  rate  (every  100  msec).  The  10  Hz 
update  of  the  cockpit  indicators  is  considered  a  hard  speci¬ 
fication,  but  some  minor  flexibility  might  be  available  for 
the  40  msec  filter  execution  time.  Assuming  40  msec  as 
worst  case  allows  a  maximum  of  12000  computer  operations  on 
the  F-4E/G's  300  KOPS  fire  control  computer  for  a  filter 
algorithm.  Another  limitation  of  this  computer  is  that  it 
has  only  8K  words  of  memory  available  for  a  Kalman  filter 
algorithm  and  supporting  auxiliary  calculations . 

The  Kalman  filter  implemented  in  this  constrained 
memory  must  satisfy  both  the  LRI  and  ACM  functions.  Unfor¬ 
tunately,  these  two  functions  operate  under  different  target 
assumptions.  LRI  assumes  a  non-maneuvering  target  and  re¬ 
quires  filter  estimates  of  target  velocity  only,  while  ACM 
assumes  a  highly  dynamic  target  and  requires  filter  esti¬ 
mates  of  target  acceleration  as  well  as  velocity  (9,10). 
Additionally,  the  noise  content  of  the  radar  measurements  is 
partially  a  function  of  target  range  and  differs  signifi¬ 
cantly  between  ACM  measurements,  taken  within  8  nautical 
miles  (nm),  and  LRI  measurements,  taken  from  8  to  30  nm  (9). 
As  a  result  of  these  different  operating  conditions,  inher¬ 
ent  differences  exist  between  the  two  functions  and  a  filter 
designed  for  LRI  might  not  prove  suitable  for  ACM.  In  the 


past,  the  Wiener-Hopf  algorithm  was  implemented  as  two  sep¬ 
arate  filters,  one  dedicated  to  ACM  and  the  other  for  LRI. 
However,  a  significant  portion  of  the  filter  computations 
are  common  to  both  functions.  By  carefully  accounting  for 
the  differences  discussed  above,  a  single  combined  filter 
algorithm  is  designed  to  run  in  either  ACM  or  LRI  mode  in 
the  latest  version  of  the  OFF.  The  combined  function  filter 
conserves  overall  memory  usage  but  exacts  a  cost  in  proces¬ 
sing  time  because  the  more  complex  combined  filter  requires 
longer  processing  time  than  either  single  function  filter 
alone.  Regardless  of  whether  separate  or  combined  filters 
are  implemented,  the  double  constraint  of  memory  space  and 
processing  time  limits  the  complexity  of  the  proposed  Kalman 
filter.  To  cope  with  this  constraint,  reduced  order  and 
constant  gain  Kalman  filters  are  investigated  in  this  study, 
A  third  constraint  imposed  by  the  computer  is  its  16-bit 
wordlength.  Short  wordlength  is  a  frequent  contributor  to 
numerical  instability  in  Kalman  filter  applications,  a  prob¬ 
lem  which  strongly  motivates  the  use  of  square  root  or  U-D 
filter  formulations  (1:368-369). 

Besides  the  constraints  imposed  by  computer  hardware, 
the  Kalman  filter  designer  must  also  consider  the  con¬ 
straints  imposed  by  the  F-4E/G  sensor  group.  As  explained 
earlier,  it  is  not  economically  feasible  at  present  to 
replace  or  upgrade  the  APQ-120  radar  or  other  onboard 
sensors  which  supply  measurements  to  the  OFF  for  fire  con¬ 
trol  purposes.  These  include  the  central  air  data  computer 


> 


(CADC)  and  inertial  navigation  system  (INS).  The  measure¬ 
ments  provided  by  these  instruments  and  the  errors  asso¬ 
ciated  with  them  must  be  taken  as  givens  for  the  purpose  of 
this  work. 

The  above  constraints  place  severe  limitations  on  a 
Kalman  filter  application  for  the  F-4E/G.  Such  a  limited 
Kalman  filter  might  not  perform  accurately  enough  to  meet 
the  needs  of  the  fire  control  system  in  a  combat  environ¬ 
ment.  To  judge  its  adequacy,  the  performance  of  the  filter 
must  be  measured  against  explicit  criteria.  In  the  case  of 
LRI,  the  filter  must  produce  target  velocity  estimates  whose 
overall  magnitude  error  has  a  maximum  standard  deviation  of 
30  feet  per  second  (8).  This  accuracy  is  required  by  the 
current  OFF  to  generate  cockpit  indications  which  meet  the 
combat  requirement  of  less  than  10%  out-of-boundary  launches 
and  less  than  10%  missed  opportunities  to  attack  (10:3-1). 

To  warrant  implementation,  the  Kalman  filter  must 
continue  to  meet  all  performance  criteria  during  any 
maneuvers  required  by  LRI  attack.  Currently,  the  Wiener- 
Hopf  filter's  calculation  of  target  velocity  is  adequate 
during  straight  and  level  flight,  but  is  virtually  useless 
during  F-4E/G  maneuvering  (9).  The  radar  azimuth  and 
elevation  positions  as  well  as  the  outputs  of  the  radar's 
rate  gyros  lag  behind  fighter  maneuvers  and  oscillate.  The 
Wiener-Hopf  filter  amplifies  this  problem  so  that  during 
maneuvers,  the  lag  and  oscillation  of  the  radar's  outputs 
become  so  severe  that  the  missile  range  estimate  collapses 
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to  zero  and  the  pilot's  steering  dot  in  the  cockpit  bounces 
erratically  across  the  screen. 


To  overcome  this  problem  and  to  minimize  pilot  distrac¬ 
tion,  a  postfilter  is  included  in  the  present  software.  The 
postfilter  essentially  puts  the  fire  control  system  into 
memory  mode  during  fighter  maneuvers  by  ignoring  the  ampli¬ 
fied  radar  noise  and  by  smoothing  the  output  to  the  steering 
dot  and  range  indicator  using  past  values  of  rate  gyro 
measurements.  This  arrangement  provides  a  non-distracting 
display  but  fails  to  provide  useful  information  to  aid  the 
pilot  in  successfully  attacking  the  target.  The  problem  is 
described  in  aircrew  technical  orders  as  follows: 

"When  the  aircraft  rolls  to  center  the  aim  dot, 
APQ-120  rate  gyro  transients  make  the  steering  dot 
jump.  To  make  the  dot  stable,  the  computer  uses 
the  rate  gyro  outputs  it  received  prior  to  roll  until 
the  fighter  roll  angle  is  established  and  the  tran¬ 
sients  settle."  (11:1-82) 

Correcting  this  problem  and  eliminating  the  need  for  post¬ 
filter  data  processing  are  the  essential  benefits  which 
could  be  derived  from  this  research.  Thus,  the  desired  end 
product  of  this  thesis  effort  is  a  Kalman  algorithm  which 
will  generate  target  velocity  estimates  of  the  accuracy 
cited  above,  regardless  of  F-4E/G  combat  maneuvering,  and 
which  can  be  implemented  in  the  onboard  computer  within  the 
cited  constraints  of  computer  memory  and  execution  time. 

Scope  and  Assumptions 

Although  integration  with  the  ACM  function  is  a  major 
consideration,  this  study  deals  directly  with  filter  design 


for  LRI  only.  Filter  algorithms  are  tested  in  a  computer 
simulation  using  a  Kalman  filter  analysis  package  called 
the  Simulation  for  Optimal  Filter  Evaluation,  or  SOFE,  and  a 
complementary  data  postprocessor  and  plotting  program  called 
SOFE  Plotter,  or  SOFEPL  (13,28).  The  Fortran  77  language  is 
used  for  the  simulation  and  no  attempt  is  made  to  implement 
the  filters  in  F-4E/G  machine  language.  This  study  does  not 
include  flight  testing,  but  recommendations  are  included  to 
assist  the  incorporation  of  the  filters  designed  into  the 
OFP  for  bench  and  flight  testing. 

Although  many  different  Kalman  filter  configurations 
are  possible  for  this  application,  only  two  filters  are 
initially  developed.  Peiformance  evaluations  are  then  used 
to  tune  these  filters  for  best  possible  performance.  If 
further  performance  improvements  are  still  required  after 
tuning,  filter  modifications  are  possible.  These  modifica¬ 
tions  include  constant  gain,  adaptive  gain,  or  the  augmenta¬ 
tion  of  additional  states  to  represent  originally  unmodelled 
effects . 

Both  of  the  initially  proposed  filter  models  are  non¬ 
linear  in  their  measurement  equations.  This  precludes  use 
of  conventional  Kalman  filter  formulations.  The  extended 
Kalman  filter  (EKF)  is  used  instead  for  these  reasons: 

(1)  its  proven  superior  performance  over  linearized 
Kalman  filters  in  cases  where  a  nominal  state 
trajectory  is  not  known  a  priori  (2:42,52  to  55, 
12:38) 


(2)  its  computational  simplicity  when  compared  to 
higher  order  estimators  (2:221,235),  and 

(3)  its  compatibility  with  the  available  performance 
evaluation  tool,  SOFE  (12:38)(13). 

Use  of  the  EKF  entails  approximations  of  model  non- 
linearities  in  some  aspects  of  the  algorithm,  and  explicit 
evaluation  of  nonlinearities  in  other  aspects,  as  explained 
in  Chapter  IV.  It  is  important  to  note  that  the  approxima¬ 
tions  inherent  in  the  EKF  formulation  might  not  be  adequate 
in  cases  where  the  nonlinearities  are  extremely  harsh 
(2:224).  The  nonlinear  models  developed  in  this  thesis  are 
assumed  to  be  sufficiently  benign  to  allow  use  of  the  EKF. 
The  validity  of  the  dynamics  models  themselves  for  repre¬ 
senting  physical  reality  is  another  implicit  assumption 
necessary  for  computer  simulation  studies.  The  models  de¬ 
veloped  in  Chapter  II  are  assumed  adequate  based  on  their 
successful  application  in  previous  tracking  problems. 

In  addition  to  the  fundamental  assumptions  above, 
several  specific  assumptions  allow  the  simplification  of 
complex  effects  which  have  relatively  minor  influence  on 
overall  system  behavior.  The  simplification  of  non-dominant 
effects  permits  this  study  to  emphasize  filter  design  and 
evaluation . 

For  LRI,  the  target  is  assumed  to  be  a  non-maneuvering 
point  mass,  while  ACM  assumes  that  targets  are  maneuvering. 
Both  modes  treat  the  F-4E/G  as  a  three  dimensional  vehicle. 


but  the  acknowledged  velocity  and  acceleration  discrepancies 
of  different  locations  on  the  airframe  during  nonzero  atti¬ 
tude  rates  are  assumed  neglible  for  LRI.  Here,  "attitude 
rate"  refers  to  any  combination  of  heading,  pitch,  or  roll 
rate.  Although  yaw  rates  have  an  admitted  effect,  they  are 
not  calculated  in  the  F-4E/G  OFF  (8).  Flight  without  side¬ 
slip  is  assumed  in  this  study.  The  OFF  calculates  fighter 
attitude  rates  by  using  the  incremental  INS  angles  of  pitch, 
roll,  and  heading.  The  angles  and  angle  rates  derived  from 
them  are  corrupted  by  sensor  and  approximation  errors. 

These  errors  are  expected  to  be  much  smaller  than  the  errors 
committed  by  the  filter,  and  are  assumed  neglible.  Although 
this  assumption  may  not  be  completely  justifiable,  the  lim¬ 
ited  memory  and  speed  of  the  fire  control  computer  preclude 
the  inclusion  of  a  complex  INS  error  model  in  the  OFF. 

The  AFQ-120  antenna  is  pointed  at  the  target  by  a 
closed  loop  control  system.  The  response  time  of  this 
system  is  assumed  to  be  very  fast  in  comparison  with  the 
filter  cycle  time.  For  the  short  time  required  to  take  a 
measurement,  the  antenna  is  assumed  to  be  space-stabilized. 
This  allows  simplification  of  the  measurement  equations 
since  the  antenna  rates  do  not  have  to  be  integrated  to 
determine  antenna  angle.  The  errors  associated  with  antenna 
angle  include  true  pointing  error  caused  by  the  control 
system  as  well  as  measurement  errors  deriving  from  resolvers 
and  other  sensors.  In  this  study,  these  errors  are  combined 


and  quantified  as  a  total  measurement  error.  Measurement 
assumptions  are  detailed  in  Chapter  III. 


Methodology 

Having  identified  the  problem,  defined  limitations,  and 
proposed  a  Kalman  filter  as  a  solution,  the  goal  of  this 
study  is  to  demonstrate  the  feasibility  of  the  proposed 
solution.  This  section  describes  the  logical  approach  used 
for  this  purpose  and  relates  each  step  to  the  remaining 
chapters  of  the  report.  Conceptually,  the  study  can  be 
broken  into  the  following  steps: 

(1)  Model  development  (Chapters  II  and  III), 

(2)  Filter  implementation  (Chapter  IV), 

(3)  Performance  evaluations  of  filters  (Chapter  V), 

(4)  Results  of  performance  evaluations  (Chapter  VI), 

(5)  Conclusions  and  recommendations  (Chapter  VII). 

These  areas  are  expanded  below. 

The  most  crucial  element  of  Kalman  filter  design  is  the 
selection  of  adequate  models  as  a  basis  for  the  filter  (1). 
In  Chapter  II,  several  dynamics  models  for  air-to-air 
tracking  are  examined  and  two  are  selected  using  the  cri¬ 
teria  of  adequacy  to  represent  reality,  computational  sim¬ 
plicity,  and  compatibility  with  the  existing  F-4E/G  fire 
control  system.  The  measurements  available  to  the  filter 
are  a  fixed  function  of  existing  F-4E/G  hardware.  The 
mathematical  modelling  of  these  measurements  is  dependent  on 


the  dynamics  models  chosen.  For  this  reason,  Chapter  III 
begins  with  an  examination  of  the  available  F-4E/G  sensor 
group  and  then  expresses  the  measurements  produced  by  these 
sensors  in  terms  of  the  state  variables  defined  for  both 
models  in  Chapter  II. 

Chapter  IV  develops  the  equations  for  an  extended 
Kalman  filter  based  on  the  models  developed  in  Chapters  II 
and  III.  This  filter  is  evaluated  using  standard  Monte 
Carlo  analysis  techniques  explained  in  Chapter  V.  That 
chapter  includes  the  development  of  a  trajectory  generation 
program  used  as  a  standard  of  comparison,  or  "truth  model", 
for  filter  performance. 

Finally,  the  results  of  the  performance  evaluations  are 
analyzed  in  Chapter  VI  and  recommendations  concerning 
further  study  or  actual  implementation  are  summarized  in 
Chapter  VII . 


1 1 .  Dynamics  Models 


Model  Selection  Process 

This  chapter  describes  practical  dynamics  models  used 
to  derive  a  Kalman  filter  for  the  F-4E/G  fire  control  sys¬ 
tem.  As  discussed  in  Chapter  I,  the  crux  of  the  Kalman 
filter  design  problem  is  obtaining  an  adequate  mathematical 
model  upon  which  to  base  the  filter  (1:289).  For  this 
reason,  considerable  time  and  thought  has  been  spent  on 
acquiring  appropriate  dynamics  models  for  this  study. 

The  general  case  of  airborne  target  tracking  by  radar 
has  been  researched  extensively  in  the  past  and  numerous 
mathematical  models  already  exist  for  this  purpose  (4,12, 
33,34,35).  These  models  typically  involve  state  variables 
defined  in  terms  of  target  position,  velocity,  and  accelera 
tion.  Table  II-l  lists  several  well  established  dynamics 
models  for  target  tracking  and  gives  references  for  each. 


TABLE  II-l 

Target  Tracking  Dynamics  Models 


Model  References 

Constant  Velocity 

1,4 

Constant  Acceleration 

4.15 

Gauss-Markov  Acceleration 

4,34 

Constant  Turn  Rate 

12,33 

Acceleration  Perpendicular 
to  Planform 

36,37 

Variations  of  the  models  listed  in  Table  ll-l  may 
include  states  defined  in  Cartesian,  spherical,  or  cylindri¬ 
cal  coordinates  and  aligned  with  inertial,  1 i ne-o f -s i gh t ,  or 
other  reference  frames.  Additionally,  the  target's  position 
velocity,  and  acceleration  states  may  be  defined  as  total 
parameters  referenced  to  a  fixed  inertial  point,  or  alter¬ 
nately,  as  relative  parameters  referenced  to  another  system 
element,  often  the  tracker.  Each  model  variation  entails 
different  assumptions,  approximations,  and  definitions. 

These  differences  can  lead  to  vastly  different  Kalman  filter 
implementations. 

For  the  specific  case  of  the  F-4E/G  fire  control  sys¬ 
tem,  the  initial  problem  is  to  select  which  variations  of 
the  above  models  are  most  appropriate  for  this  particular 
application.  In  this  study,  the  model  selection  process 
addresses  three  issues: 

1.  adequacy  of  the  model  to  describe  the  physical 
behavior  of  the  system, 

2.  complexity  of  the  model,  and 

3.  compatibility  of  the  model  with  existing  F-4E/G 
fire  control  system  hardware  and  software. 

The  first  consideration  is  the  adequacy  of  the  proposed 
models  to  describe  the  actual  physical  behavior  of  antici¬ 
pated  airborne  targets.  In  reality,  the  F-4E/G's  targets 
are  various  types  of  military  aircraft.  Such  targets  under¬ 
go  unpredictable  accelerations  with  resulting  changes  in 


velocity  and  position.  The  accelerations,  although  unpre¬ 
dictable  in  both  direction  and  magnitude,  are  bounded  by  the 
structural  g-limitations  of  the  particular  aircraft  type. 
Small,  highly  maneuverable  fighter  aircraft  are  capable  of 
higher  frequency  variations  of  acceleration  over  a  larger 
range  of  extremes  than  large,  less  maneuverable  bomber  or 
transport  aircraft.  Either  type  of  target  is  well  described 
by  band-limited  random  acceleration  modelled  as  a  first 
order  Gauss-Markov  process  (6:157): 

da  1 

--  =  -  a  +  w  (2-1) 

dt  T 

where  ^  is  the  target  acceleration  vector,  i  is  a  time- 
constant  of  target  acceleration  change  associated  with  the 
type  of  aircraft,  and  w  is  a  zero-mean,  white  Gaussian  noise 
vector.  Other  models  besides  the  Gauss-Markov  model  can 
accurately  describe  the  motion  of  randomly  accelerating 
aircraft  targets.  Constant  turn  rate  and  planf orm-dependent 
acceleration  models  are  examples  (12,33,36).  Conversely, 
constant  velocity  or  constant  acceleration  models  are  not 
usually  good  models  for  randomly  accelerating  targets 
because  they  assume  a  target  with  either  zero  acceleration 
or  a  steady  maneuver.  However,  in  the  case  of  the  F-4E/G 
LRI  function,  the  effect  of  extended  target  range  is  that 
random  accelerations  caused  by  target  maneuvers  are  obscured 
by  system  noise.  Indeed,  the  current  LRI  function  assumes  a 
non-maneuvering  target  and  does  not  use  filter-computed 
estimates  of  target  acceleration  (10).  For  this  reason, 


constant  velocity  or  acceleration  models  might  be  adequate 
for  the  LRI  function,  but  are  not  acceptable  for  ACM. 

A  second  consideration  in  selecting  an  appropriate 
dynamics  model  for  the  F-4E/G  filter  is  model  complexity. 
Complexity  becomes  an  issue  when  hardware  constraints  limit 
the  amount  of  computer  memory  available  for  filter  implemen¬ 
tation  and  computer  speed  restricts  the  amount  of  computa¬ 
tion  possible  during  real  time  operation.  As  discussed  in 
Chapter  I,  these  constraints  are  severe  in  the  F-4E/G  fire 
control  system  computer.  The  possible  8K  words  of  memory 
available  for  the  LRI  filter  shrink  significantly  when 
space  is  reserved  for  complex  coordinate  transformations, 
residual  monitoring  algorithms,  adaptive  estimation  tech¬ 
niques,  or  decision  logic  for  activating  either  ACM  or  LRI 
filters.  Using  Fig.  2-1  as  a  guide,  the  maximum  number  of 
states  possible  for  a  Kalman  filter  implementation  in  the 
available  F-4E/G  computer  memory  is  estimated  to  be  between 
12  and  15.  The  combination  of  slow  speed  (300  KOPS)  and 
limited  time  available  for  each  filter  recursion  imposes 
an  even  more  severe  constraint  and  restricts  the  number  of 
filter  states  to  the  lower  end  of  this  range.  In  final 
analysis,  the  severe  constraints  of  the  F-4E/G  fire  control 
system  hardware  dictate  the  choice  of  the  simplest  filter 
possible  which  meets  performance  specifications.  This 
situation  is  not  uncommon  in  engineering  applications  and  is 
often  espoused  as  a  general  design  philosophy  in  textbooks 
(1:290,356). 


(1:355) 

Fig.  2.1,  Kalman  Filter  Computer  Memory  Requirements 

In  general,  the  complexity  of  a  given  Kalman  filter 
implementation  is  both  model  dependent  and  geometry  dependent 
The  choice  of  a  particular  dynamics  model  often  determines 
the  number  of  filter  states  required.  For  instance,  a 
constant  velocity  model  typically  uses  fewer  states 
than  a  model  which  includes  several  acceleration  states. 


Although  the  number  of  states  is  an  important  determinant, 
filter  complexity  is  also  a  fraction  of  model  nonlinearities 


and  required  auxiliary  computations  such  as  coordinate 
transformations.  These  factors  are  closely  related  to  the 
coordinate  reference  used  to  define  the  filter  states.  For 
example,  a  Gauss-Markov  acceleration  model  results  in  linear 
dynamics  equations  when  implemented  in  inertial  Cartesian 
coordinates.  The  same  model  defined  in  a  Cartesian  line-of- 
sight  reference  has  nonlinear  dynamics  (15,16).  The  increased 
complexity  and  computational  burden  resulting  from  model 
nonlinearity  is  discussed  in  greater  detail  in  Chapter  IV. 

Given  the  discussion  to  this  point,  a  constant  velocity 
model  defined  in  a  reference  frame  which  yields  simple 
linear  dynamics  would  be  a  good  choice  for  the  F-4E/G's  LRI 
filter  because  it  is  the  simplest  model  of  those  postulated 
as  adequate.  However,  a  third  issue  must  be  addressed 
before  the  final  choice  is  made.  The  LRI  function  does  not 
operate  in  isolation,  but  in  conjunction  with  the  total  OFF. 
The  LRI  filter  must  be  compatible  with  the  ACM  function  and 
other  OFF  components.  Also,  it  must  make  use  of  the  exist¬ 
ing  measurements  from  the  available  F-4E/G  sensor  group,  as 
will  be  discussed  in  Chapter  III.  The  ACM  function  does  not 
share  the  LRI  assumption  of  a  non-maneuvering  target.  In¬ 
stead,  it  assumes  a  highly  dynamic  target  at  close  range  and 
requires  filter  computed  estimates  of  the  target's  accelera¬ 
tion  (10).  Thus,  the  ACM  filter  demands  a  more  complex 
dynamics  model  than  LRI  alone  and  any  combined  ACM/LRI 
filter  implementation  for  the  F-4E/G  OFF  must  take  this  into 


account . 


Using  the  criteria  of  adequacy,  simplicity,  and  compat¬ 
ibility  as  discussed  above,  two  candidate  filters  are 
initially  proposed,  one  based  on  a  constant  velocity  (CV) 
model  and  the  other  based  on  a  Gauss-Markov  acceleration 
(GMA)  model.  These  models  are  developed  in  the  following 
sections . 

Constant  Velocity  ( CV )  Model 

A  six-state  constant  velocity  model  defined  in  Carte¬ 
sian  inertial  coordinates  is  proposed  as  the  basis  of  an 
acceptable  LRI  filter.  This  filter  represents  a  solution  of 
minimal  complexity  and  is  chosen  based  on  the  philosophy  of 
using  the  simplest  possible  filter  that  yields  adequate 
results  (1:290,  6:159,  17:423).  In  addition  to  minimal 
state  dimension,  this  model  has  the  advantage  of  linear, 
time-invariant  dynamics  equations.  These  equations  lead  to  a 
constant-valued  state  transition  matrix  in  the  filter  imple¬ 
mentation,  assuming  the  filter  measurement  interval  is  con¬ 
stant.  This  greatly  simplifies  filter  propagation,  as 
explained  in  Chapter  IV, 

The  major  drawback  of  the  CV  model  is  that  it  does  not 
model  the  target  acceleration,  which  is  a  requirement  for 
ACM.  In  the  F-4E/G  OFF,  a  CV-based  filter  would  either 
be  implemented  as  a  separate  LRI  filter  or  integrated 
into  a  more  complex  ACM  filter.  In  the  context  of  this 
study,  the  CV  filter  serves  as  an  initial  demonstration  of 
feasibility,  as  well  as  a  viable  LRI-only  filter. 


Fig.  2.2.  Geographic  Inertial  Reference  Frame 

The  state  variables  for  the  CV  model  are  defined  as  the 
relative  position  of  the  target  with  respect  to  the  F-4E/G 
tracker  and  the  total  velocity  of  the  target  with  respect  to 
an  inertial  reference.  For  this  study,  the  inertial  refer¬ 
ence  is  a  geographic  frame  fixed  to  the  earth's  surface. 

The  three  Cartesian  reference  axes  for  this  frame  are  ori¬ 
ented  to  north,  east,  and  down  as  depicted  in  Fig.  2.2.  The 
effect  of  earth  rotation  is  assumed  neglible  for  the  short 
duration  encounters  (1-2  minutes)  and  comparatively  short 
distances  (8-20  miles)  anticipated  in  the  LRI  combat  environment 


Two  considerations  motivate  the  choice  of  relative 


target  position  states.  First,  relative  position  states  are 
limited  in  magnitude  to  the  value  of  the  target's  actual 
range  from  the  tracker,  whereas  total  position  states  can 
exhibit  much  larger  magnitudes  in  an  earth  fixed  reference 
frame.  Table  II-2  contains  typical  values  for  comparison. 

TABLE  II-2 


Maximum  Magnitudes  of  Relative  vs.  Total  States 


State 

Relative 

Total 

Position  (feet) 

182,400® 

4,152,640^ 

Velocity  (fps) 

1520^^ 

760^^ 

^  based  on  30  nm  required  LRI  tracking  range  (9) 
based  on  683  nm  combat  radius  (20:375) 

Q 

based  on  nominal  760  feet  per  second  (fps)  combat 
airspeed  (10:3-417) 

Numerical  stability  of  the  derived  filter  is  improved  by 
keeping  the  magnitudes  of  the  position  states  as  small  as 
possible,  since  small  differences  of  large  numbers  are  a 
common  source  of  numerical  instability  in  Kalman  filters 
(1:256).  Keeping  the  magnitudes  of  the  states  small  reduces 
the  likelihood  of  this  type  of  instability.  Excessive  re¬ 
duction  of  state  magnitudes,  though,  may  result  in  filter 
outputs  which  lack  sufficient  precision.  Such  is  the  case 
when  very  large  units  are  used  in  an  attempt  to  reduce  the 


position  state  magnitudes.  Target  position  estimates  ex¬ 
pressed  in  miles,  for  instance,  are  insufficiently  precise 
for  F-4E/G  fire  control  system  functions.  For  this  reason, 
the  OFF  uses  feet  as  the  standard  unit  of  distance.  Expres 
sing  position  in  feet  requires  14  bits  to  represent  the 
largest  relative  position  state  magnitude  in  Table  II-2  to 
within  eight  foot  precision.  It  requires  18  bits  to  repre¬ 
sent  the  largest  total  position  state  to  the  same  precision 
(see  Appendix  A).  Consequently,  relative  states  produce 
position  magnitudes  which  achieve  the  best  compromise  be¬ 
tween  numerical  stability  and  precision.  Note  that  in  a 
computer  simulation,  the  magnitude  of  total  position  states 
could  be  sufficiently  limited  by  initializing  the  attacker 
at  an  inertial  origin  each  time  and  limiting  the  simulation 
time  to  very  short  periods.  However,  this  is  not  practical 
for  the  real  world  environment  and  thus,  relative  position 
states  are  again  desirable. 

A  second,  more  important  reason  to  use  relative  posi¬ 
tion  states  is  the  nature  of  the  available  F-4E/G  target 
measurements.  Radar  measurements  are  inherently  relative 
and  require  the  incorporation  of  attacker  ownship  inertial 
position  in  order  to  generate  measurements  of  the  target's 
inertial  position.  Using  relative  target  position  states 
eliminates  the  need  for  F-4E/G  ownship  position.  Since 
ownship  position  information  is  not  available  to  the  OFF  in 
the  F-4E/G  fire  control  system,  relative  position  states 
are,  in  fact,  necessary. 


Although  similar  considerations  affect  the  choice  of 


velocity  states,  the  driving  factor  here  is  the  enhancement 
of  the  modelling  assumptions.  As  with  position  states,  it 
is  desirable  from  a  numerics  standpoint  to  limit  the  magni¬ 
tudes  which  the  velocity  states  can  realize.  Neglecting 
wind,  inertial  target  velocity  is  limited  by  the  target's 
maximum  airspeed  capability.  Using  relative  velocities 
essentially  doubles  this  range  of  values,  with  extremes 
occurring  in  cases  of  directly  opposed  trajectories. 

Table  II-2  gives  typical  values  for  relative  and  total 
velocities.  In  this  case,  7  bits  achieve  a  4  fps  precision 
for  the  maximum  total  velocity  value  of  760  fps,  while  8 
bits  achieve  similar  precision  for  the  maximum  relative 
velocity  of  1520  fps  (see  Appendix  A).  Since  it  requires 
only  one  additional  bit  to  obtain  similar  precision,  the 
effect  of  using  either  of  these  velocity  ranges  is  essen¬ 
tially  the  same . 

A  more  critical  consideration  is  again  compatibility 
with  existing  OFP  software.  The  F-4E/G's  LRI  function  is 
presently  designed  to  use  total  target  velocity  estimates 
from  the  filter.  This  fact  in  itself  provides  a  strong 
incentive  to  use  total  velocity  states.  Further,  in  the  CV 
model,  the  velocity  states  are  assumed  constant.  The  LRI 
function  assumes  a  non-manuever ing  target,  which  supports 
the  model  of  constant  target  inertial  velocity.  However, 
the  F-4E/G  itself  must  be  assumed  to  be  maneuvering  into  an 
optimal  attack  position.  Since  the  attacker  maneuvers,  the 
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relative  velocity  is  not  constant  and  the  validity  of  the  CV 
model  is  degraded  by  using  relative  target  velocity  states. 

Clearly,  the  use  of  relative  position  and  total  veloc¬ 
ity  states  is  viable,  but  this  combination  does  give  rise  to 
some  complicating  factors.  One  consequence  of  using  total 
target  velocity  states  in  the  same  model  as  relative  posi¬ 
tion  states  is  that  the  F-4E/G’s  ownship  velocity  must  be 
calculated  and  subtracted  from  the  position  state  dynamics 
equations  and  certain  measurement  equations.  Fortunately, 
this  does  not  create  a  problem  in  the  F-4E/G  because  ownship 
velocity  is  already  calculated  in  the  current  OFF  from  true 
airspeed  and  angle-of-attack  inputs. 

Using  the  states  developed  above,  the  CV  model  state 
vector  is : 


where  x^(t),  x^(t),  and  x^(t)  in  (2-2)  are  defined  as  the 

components  of  target  position  relative  to  the  F-4E/G  tracker 
in  the  north,  east,  and  down  directions  respectively, 

(  labeled  x^,  x^.,  and  x^  in  Fig.  2.2);  and  x^,  x^,  and  x^  are 
the  respective  north,  east,  and  down  components  of  the 
target's  total  inertial  velocity. 


The  following  set  of  differential  equations  governs  the 


dynamics  of  the  CV  state  vector: 

x^(t)  =  x^Ct)  -  vjj^(t)  (2-3) 

X2(t)  =  0  +  Wj^(t)  (  2-4) 

X3(t)  =  x^(t)  -  vjj,(t)  (2-5) 

x^( t )  =  0  +  Wg( t )  ( 2-6) 

X3(t)  =  x^(t)  -  V^j^(t)  (2-7) 

x^(t)  =  0  +  (2-8) 


where  ^p(t)  is  a  vector  containing  F-4E/G  ownship  inertial 
velocity  at  time  t.  Superscript  I  and  subscript  F  denote 
inertial  F-4E/G  quantities,  while  the  subscripts  N,  E,  and  D 
indicate  north,  east,  and  down  components,  respectively. 

The  vector  w ( t )  contains  the  north,  east,  and  down  compo¬ 
nents  of  zero-mean,  white,  Gaussian  driving  noise  which  are 
included  Co  account  for  inherent  modelling  errors. 

Errors  exist  in  all  models,  especially  highly  simpli¬ 
fied  versions  such  as  the  CV  model.  These  errors  exist 
because  all  factors  which  affect  the  dynamics  of  a  real 
world  system  are  not  included  in  the  simple  model.  A  re¬ 
duced  order  model  typically  includes  only  the  dominant 
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factors  which  affect  the  system.  To  account  for  the  uncer¬ 
tainty  introduced  by  modelling  errors,  the  components  of  the 
driving  noise  vector,  w ( t ) ,  are  added  to  the  model  equa¬ 
tions.  The  strength  of  this  noise,  Q(t),  is  a  function  of 
the  various  approximations  and  unmodelled  effects  associated 
with  a  given  model. 

In  the  CV  model,  target  velocity  is  assumed  constant, 
yet  all  real  targets  undergo  some  acceleration  and  velocity 
fluctuations.  Even  an  aircraft  attempting  to  fly  straight 
and  level  at  constant  airspeed  experiences  minor  velocity 
variations  caused  by  complex  aircraft  dynamics  and  a  variety 
of  atmospheric  conditions.  At  distant  ranges,  these  small 
fluctuations  are  virtually  undetectable  and  suggest  that  the 
terms  of  Q(t)  are  small.  For  closer  ranges  or  a  maneuvering 
target,  Che  terms  of  Q(t)  might  be  changed  to  account  for 
greater  uncertainty  in  the  model. 

Another  unmodelled  effect  is  earth  rotation.  Assumed 
neglible  for  brief  encounters  at  short  ranges,  earth  rota¬ 
tion  becomes  an  increasingly  significant  source  of  error  as 
target  range  and  Cracking  time  increase.  Accordingly,  Che 
value  of  Q(t)  due  to  this  effect  increases  with  range  and 
tends  to  offset  the  error  described  above. 

A  third  source  of  error  in  the  CV  model  is  the  error 
associated  with  F-4E/G  ownship  velocity.  Mathematically, 
the  velocity  required  in  (2-3),  (2-5),  and  (2-7)  is 

r 

the  inertial  velocity  of  the  tracking  radar  in  Che  nose  of 
Che  aircraft.  Realistically,  the  Xp(t)  supplied  to  the  OFF 


IS  an  approximation  based  on  noise-corrupted  measurements  of 
true  airspeed  and  angle-o f -at tack .  There  are  measurement 
errors  associated  with  the  instruments  which  measure  true 
airspeed  and  angle-of-attack ,  but  within  the  limits  of  these 
errors,  the  approximation  for  F-4E/G  inertial  velocity  is 
valid  for  any  point  on  the  aircraft,  providing  the  F-4E/G 
has  no  attitude  rate.  When  an  attitude  rate  is  present,  the 
approximation  is  still  valid  for  points  lying  along  the 
aircraft's  instantaneous  axis  of  rotation,  but  becomes  in¬ 
creasingly  erroneous  at  points  further  from  this  axis. 

Since  the  F-4E/G  tracking  radar  is  located  well  forward  of 
the  aircraft's  usual  attitude  rotation  axis,  it  is  reason¬ 
able  to  assign  increased  error  to  during  nonzero  F-4E/G 
attitude  rates  and  to  adjust  Q(t)  accordingly.  Correction 
factors  for  the  errors  in  caused  by  F-4E/G  attitude 

rates  add  significant  complexity  to  the  dynamics  model. 

These  factors  are  not  included  in  this  study  because  the 
constraints  of  the  F-4E/G  fire  control  computer  require 
minimal  complexity. 

The  above  errors  are  combined  into  a  single  vector 
stochastic  process,  w ( t  )  .  By  virtue  of  the  central  limit 
theorem,  t )  is  a  zero-mean,  white,  Gaussian  process, 
regardless  of  the  probability  distributions  of  the  component 
errors  (1:109).  Q(t)  is  then  eX;ressed  by  the  covariance 
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Actual  values  for  Q(t)  are  difficult  to  determine  a  priori. 
For  this  study,  a  nominal  value  of  5000  ft  /sec  is  chosen 
as  a  reasonable  first  guess  for  the  diagonal  elements  of 
Q(t).  This  value  can  be  increased  or  decreased  to  improve 
filter  performance,  as  discussed  in  Chapter  V.  Several 
weaknesses  in  this  noise  model  exist.  First,  noise 
strengths  are  assumed  constant  in  this  model.  This  is  not 
true  in  reality  because  resolving  target  velocity  into 
north,  east,  and  down  components  yields  noise  strengths 
which  are  dependent  on  changing  angular  rates  and  accelera¬ 
tions.  Additionally,  the  components  of  the  driving  noise 
are  modelled  as  independent  of  each  other  and  the  off- 
diagonal  terms  of  Q(t)  are  assumed  zero.  This  aspect  of  the 
noise  model  is  questionable  because  some  degree  of  correla¬ 
tion  among  the  noise  components  is  unavoidable.  These 
issues  are  not  resolved  in  this  study  and  provide  an  area 
for  further  investigation  if  the  performance  of  the  basic  CV 
model  does  not  meet  the  criteria  outlined  in  Chapter  I  for 
the  LRI  function. 

If  the  LRI  function  were  the  only  function  of  the  OFF 
requiring  filtered  estimates  of  target  parameters,  the  CV 
model  ji(?ne  might  represent  an  adequate  filter  using  minimal 
memory  and  processing  time.  However,  as  discussed  earlier, 
the  F-4E/G  OFF  also  has  a  requirement  for  target  estimation 
in  Che  ACM  environment  and  this  needs  a  more  complex  model 
than  the  simple  CV  model.  For  this  reason,  a  Gauss-Markov 
acceleration  model  is  also  proposed. 


A  nine-state  Gauss-Markov  acceleration  model  defined  in 
Cartesian  inertial  coordinates  is  proposed  as  an  adequate 
filter  model  for  both  ACM  and  LRI  functions.  This  GMA  model 
represents  a  modest  increase  in  complexity  when  compared  to 
the  CV  model,  but  still  has  an  acceptably  low  state  dimen¬ 
sion  for  OFF  implementation  (see  Fig.  2,1).  The  state 
variables  chosen  for  the  GMA  model  are  total  target  accel¬ 
eration,  total  target  velocity,  and  target  position  relative 
to  the  F-4E/G  tracker.  All  states  are  coordinatized  in  the 
geographic  inertial  coordinate  system  depicted  in  Fig.  2,2, 
which  is  identical  to  the  frame  used  in  the  CV  model.  As  in 
the  CV  model,  the  use  of  this  non-rotating  coordinate  frame 
leads  to  linear,  time-invariant  dynamics  equations  and  the 
resultant  filter  has  benefits  as  described  earlier,  provided 
the  time  step  between  filter  iterations  is  constant. 

The  GMA  state  vector  is  expressed  as: 


where  x^(t),  x^(t),  and  x^(t)  are  the  respective  north, 
east,  and  down  components  of  target  position  relative  to  the 
F-4E/G  tracker  x^,  and  Xp  in  Fig.  2.2);  X2(t),  x^(t), 

and  Xg(t)  are  north,  east,  and  down  components  of  target 
inertial  velocity;  and  x^(t),  Xg(t),  and  x^Ct)  are  north, 
east,  and  down  components  of  target  inertial  acceleration. 

The  GMA  model  is  based  on  the  assumption  that  the  time 
history  of  total  target  acceleration  is  well  modelled  by  the 
first  order  Gauss-Markov  process  expressed  in  (2-1).  This 
model  has  been  used  effectively  for  target  tracking  problems 
similar  to  the  F-4E/G  fire  control  problem  (4,  6,  12). 
Previous  work  also  indicates  that  this  model  is  better 
suited  for  modelling  total  target  acceleration  than  for 
describing  the  relative  acceleration  between  a  target  and 
airborne  tracker  (15,  18,  19).  For  this  reason,  total 
target  acceleration  states  are  used  for  the  GMA  filter. 

Total  velocity  states  are  also  a  logical  choice  for  the 
GMA  filter.  As  discussed  in  the  CV  model  section,  total 
velocity  states  enhance  filter  numerics  slightly  and  also 
have  the  advantage  of  being  in  a  form  which  is  directly 
usable  in  other  OFP  functions.  The  enhancement  of  the  CV 
modelling  assumption  is  not  a  factor  here,  but  total  veloc¬ 
ity  is  consistent  with  the  GMA  model  of  total  acceleration. 
Using  total  velocity  states  with  total  acceleration  states 
eliminates  the  need  to  correct  for  ownship  acceleration  in 
the  GMA  differential  velocity  equations,  (2-12),  (2-15),  and 
(2-18)  below.  This  is  advantageous  since  forthcoming 


versions  of  the  OFP  will  phase  out  calculations  of  ownship 
acceleration  (8).  Finally,  total  velocity  states  facilitate 
the  comparison  of  CV  and  GMA  model  performance  for  the 
purpose  of  this  study. 

Relative  position  states  are  used  in  the  GMA  model  for 
the  identical  reasons  that  relative  states  are  used  in  the 
CV  model.  As  discussed  earlier,  relative  position  states 
enhance  numerical  stability  and  are  compatible  with  avail¬ 
able  radar  measurements  and  the  available  F-4E/G  ownship 
velocity  terms  needed  as  correction  factors  in  the  differ¬ 
ential  position  equations,  (2-11),  (2-14),  and  (2-17)  below. 

For  the  GMA  dynamics  model,  the  differential  equations 
which  follow  from  the  above  definitions  are: 


x^(t)  =  X2(t)  -  Vp^(t) 

x^(t)  =  x^( t  ) 

1 

Xj(t)  =  -  -  X^(t)  +  Wj^(t) 

x^(t)  =  x^(t)  -  Vp^(t) 

X  ^  (  t )  =  X  ^  (  t  ) 


(2-11) 

(2-12) 

(2-13) 

(2-14) 

(2-15) 


1 

=  -  7  x.(t)  +  w  (t) 


x,(t) 


(2-16) 


(2-17) 


Xy(t)  =  Xg(t)  -  Vpjj(t) 

Xg(t)  =  Xg(t)  (2-18) 

1 

”  "r  (2-19) 


where  ^p(t)  is  vector  ownship  velocity,  ^( t )  is  a  vector  of 
zero-mean,  white,  Gaussian  driving  noise,  and  T  is  a  time 
constant  associated  with  target  acceleration  capability. 

The  velocity  V.p(*^)  identical  to  the  ownship  velocity  in 
the  CV  model.  The  strength  of  the  driving  noise  is  depen¬ 
dent  on  the  modelling  errors  peculiar  to  the  GMA  model. 

Some  errors  such  as  the  effect  of  earth  rotation  and  atti¬ 
tude  rate  induced  errors  associated  with  ^p(t)  are  the  same 
as  in  the  CV  model.  However,  structural  differences  in  the 
two  models  lead  to  different  strengths  of  driving  noise. 

The  GMA  model  accounts  for  target  accelerations  and  thus 
allows  changing  target  velocity.  How  close  the  GMA  model 
characterizes  target  acceleration  capability  is  determined 
largely  by  the  accuracy  of  the  time  constant  of  target 
maneuver,  t.  The  value  of  x  is  a  function  of  the  type  of 
target  aircraft  and  of  changing  flight  conditions  such  as 
true  airspeed  and  altitude.  Since  the  value  of  i  varies 
with  different  targets  and  different  conditions,  the  selec¬ 
tion  of  any  fixed  value  for  x  leads  to  some  degree  of  error. 
For  aircraft,  x  falls  into  a  fairly  narrow  range,  though. 


and  can  be  approximated  accurately  for  various  classes  of 
targets.  Highly  maneuverable  military  fighter  aircraft 
typically  exhibit  t  values  on  the  order  of  0.5  seconds  (15), 
which  is  the  value  used  in  this  study.  Less  maneuverable 
transport  aircraft  may  have  x  values  of  2-3  seconds  (15). 

For  some  applications,  the  values  of  x  associated  with  the 
different  axes  of  a  Cartesian  coordinate  system  may  not  be 
identical.  In  this  study,  however,  x  is  given  the  same 
value  for  each  direction  because  the  target  is  assumed 
capable  of  any  flight  attitude  and  therefore  equally  capable 
of  acceleration  in  any  direction. 

The  value  of  x  is  an  important  consideration  in  deter¬ 
mining  a  reasonable  value  of  Q(t)  for  the  GMA  model.  As  in 
the  CV  model,  however,  initial  estimates  of  Q(t)  are  diffi¬ 
cult  to  evaluate  with  certainty.  A  nominal  value  of  5000 
ft‘*/sec®  is  used  as  the  initial  value  of  Q(t)  in  this  study. 
This  value  is  not  directly  comparable  to  the  value  chosen 
for  the  CV  model  because  ^( t )  represents  different  effects 
in  the  two  models.  Additionally,  the  units  are  different 
from  those  used  for  noise  entering  at  the  velocity  level,  as 
in  the  CV  model.  Both  values  of  Q(t)  are  considered  nominal 
starting  points  for  the  filter  tuning  procedure  discussed 
later  . 

The  GMA-based  filter  is  expected  to  produce  target 
velocity  and  acceleration  estimates  which  meet  the  perform¬ 
ance  criteria  detailed  in  Chapter  I.  If  its  actual  perform¬ 
ance  does  not  meet  this  criteria  in  the  computer  simulation. 


other  models  might  be  considered,  including  higher  order 
models  . 


Summary 

This  chapter  has  defined  dynamics  models  for  a  Kalman 
filter  suitable  for  the  F-4E/G  fire  control  system.  Each 
model  described  has  well  defined  state  variables.  In  the 
next  chapter,  the  available  measurements  on  the  F-4E/G  are 
examined  and  mathematically  modelled  in  terms  of  the  state 
variables  just  defined. 
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III.  Measurement  Models 

Measurement  Feedback  in  Kalman  Filters 

Real  world  estimators  based  on  simplified  mathematical 
models  require  feedback  to  control  uncertainties  associated 
with  those  models.  Each  dynamics  model  developed  in  the 
previous  chapter  potentially  provides  a  basis  for  an  open 
loop  estimator  of  F-4E/G  target  parameters.  Such  an  estima¬ 
tor  would  rapidly  accumulate  large  errors  caused  by  random 
inputs  from  the  environment  and  by  inherent  deficiencies  of 
the  model  itself.  To  prevent  errors  from  growing  unbounded, 
practical  estimators  employ  feedback  of  measured  quantities 
from  the  real  dynamic  system  (7:421).  In  Kalman  filter 
algorithms,  the  feedback  of  measurements  is  implemented  in  a 
predictor-corrector  structure.  Available  measurements  are 
modelled  mathematically  in  terms  of  state  variables  from  the 
dynamics  equations.  This  model  is  included  in  the  filter 
algorithm  and  produces  a  prediction  of  the  measurement  vec¬ 
tor  for  each  measurement  time.  The  predicted  measurements 
are  subtracted  from  the  actual  measurements  to  form  a  resid¬ 
ual  (innovation)  which  is  in  turn  multiplied  by  an  optimal 
weighting  matrix  to  form  a  correction  factor.  This  factor 
is  applied  to  the  current  output  of  the  dynamics  model  to 
produce  a  corrected  state  estimate  (1:217,218). 

Deriving  an  accurate  measurement  model  requires  extensive 
knowledge  of  the  physical  nature  of  available  measurements. 


In  the  F-4E/G,  the  available  measurements  are  limited  to 
those  produced  by  existing  equipment  in  the  fire  control 
system.  Consequently,  this  chapter's  discussion  of  measure¬ 
ment  models  begins  with  a  survey  of  the  hardware  comprising 
the  available  F-4E/G  sensor  group.  The  measurements  pro¬ 
duced  by  this  set  of  existing  sensors  are  analyzed  and 
catalogued.  The  section  following  this  hardware  survey  then 
develops  the  actual  equations  used  to  model  each  measurement 
in  terms  of  CV  or  GMA  state  variables.  Finally,  the  chapter 
concludes  with  a  characterization  of  the  errors  associated 
with  each  measurement. 

Available  F-4E/G  Sensor  Group 
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antenna  is  gimbal  mounted,  allowing  ±  60°  of  travel  in  both 
azimuth  and  elevation.  The  radar  has  two  pulse  repetition 
rates  (prf),  one  for  long  range,  low  resolution  operation, 
and  the  other  for  short  range,  high  resolution  operation. 
These  and  other  detailed  performance  specifications  are 
given  in  Reference  21,  Section  2. 

The  radar  subsystem  consists  of  over  20  separate 
components,  designated  as  line  replaceable  units  (LRU's). 

A  full  listing  of  radar  LRU's  is  available  in  Reference  21, 
Section  5.  Table  III-l  lists  important  LRU's  which  produce 
or  process  radar  measurements. 

TABLE  III-l 

APQ-120  Radar  Line  Replaceable  Units 


LRU  Number 

Description 

1 

4 

7 

16 

17 

• 

Digital  ACM  computer 

Gyro-stabilized  platform 

Antenna  control  (servo  assembly) 

Antenna,  including  gimbal  assembly,  feedhorn, 
reflector,  and  rate  gyros 

Electrical  synchronizer 

Although  not  strictly  a  radar  component,  the  fire  control 
computer  performs  radar  control  functions  in  some  modes  and 
is  included  as  a  radar  LRU.  These  control  functions  are  not 
related  to  the  filter  function  of  the  computer  (20,  21). 

Various  LRU's  produce  different  measurements,  a  fact  that 
has  significant  implications  for  error  models.  For  example, 
a  single  circuit  in  the  electrical  synchronizer  (LRU  17) 


produces  both  the  range  (R)  and  range  rate  (fi)  measurements 
(21:2-19,  2-143).  This  leads  to  correlated  noise  corruption 
between  R  and  R  measurements  as  discussed  subsequently  in 
this  chapter.  Conversely,  separate  sets  of  induction  poten¬ 
tiometers  and  resolvers  in  the  antenna  gimbal  assembly  pro¬ 
duce  the  antenna  azimuth  and  elevation  angle  signals,  X  and 
.  Also,  separate  rate  gyros  measure  the  antenna  spatial 

rate  components,  to  ,  azimuth  rate,  and  to  ,  elevation  rate 

^  J 

(21:2-19,3-341).  Although  common  power  plant  and  mechanical 

interconnection  contribute  some  noise  correlation,  these 

devices  are  assumed  sufficiently  separated  to  produce  four 

measurements  which  are  independent  of  each  other  and  of  the 

range  and  range  rate  signals  as  well. 

The  F-4E/G  INS  is  primarily  a  navigation  device,  but  also 

provides  measurements  of  aircraft  pitch  (9),  roll  ( tt*)  ,  and 

heading  {'V)  angles  to  the  fire  control  system.  These  angles 

do  not  directly  provide  target  information,  but  allow  radar 

measurements  made  in  radar  reference  (i  ,i  ,k  )  coord- 

o  o  o 

inates  to  be  transformed  to  geographic  (N,E,D)  coordinates 
(see  Appendix  C).  The  CADC  uses  inputs  from  the  F-4E/G 
pitot-static  system  and  a n g 1 e-o f -a t t ac k  sensors  to  calculate 
altitude,  true  airs.eed  (TAS),  and  stall  warning  signals  for 
the  cockpit.  The  CADC  TAS  and  true  a  ng  1  e-o  f -a  1 1  ac  k  (  Qi) 
measurements  also  provice  the  OFF  with  the  ownship  vel- 
M  itv  aiding  terms  required  in  the  Kalman  filter  dynamics 
models.  Equations  (2-3),  (2-5),  (2-7)  and  (2-11),  (2-14), 
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In  addition  to  the  measurements  above,  the  radar,  INS, 
CADC,  and  other  aircraft  systems  send  and  receive  numerous 
measurements  and  control  signals  to  the  fire  control  com¬ 
puter.  Fig.  3.2  shows  the  overall  flow  of  information  in 
the  fire  control  system.  Most  of  these  signals  are  of 
peripheral  interest  in  this  study  and  are  included  here  only 
to  provide  context.  The  focus  remains  on  the  measurements 
already  described. 

The  format  in  which  Che  computer  processes  and  stores 
these  measurements  is  important.  Prior  to  storage,  the 
analog  sensor  signals  are  conditioned  by  buffering,  prescal¬ 
ing,  multiplexing,  and  A/D  conversion.  Each  of  these  opera¬ 
tions  is  an  additional  source  of  measurement  noise  corrup¬ 
tion,  as  explained  later.  The  conditioned  measurements 
leave  the  A/D  converter  appropriately  scaled  and  are  loaded 
directly  into  core  memory  (20:2-50,2-66).  Some  measurements 
are  giveii  two  different  scalings  and  stored  as  two  separate 
values.  Range,  for  instance,  is  stored  in  one  location  as 
Coarse  Range  with  a  least  significant  bit  (LSB)  value  of  232 
feet  per  bit  (air-to-air  mode),  and  stored  in  another  loca¬ 
tion  as  Fine  Range  with  a  LSB  valued  at  58.4  feet  per  bit 
(20:2-279  ).  All  angle  measurements  (  ,  9,  il>  ,  4*)  are 

also  stored  as  two  separate  measurements,  one  representing 
the  sine  of  the  angle  and  the  other  the  cosine.  Table  III-2 
lists  all  measurements  available  to  the  fire  control  compu¬ 
ter  as  they  are  actually  scored  in  Che  computer  memory. 
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TABLE  III-2  (sheet  1  of  2) 
Measurement  Storage  in  Computer  Memory 


DMA 

LOCATION 

(OCTAL) 


INPUT  SOURCE 


Clock 

Discrete  inputs 

Discrete  inputs 

Spare 

Spare 

Spare 

Spare 

Spare 

Spare 

Spare 

Spare 

Spare 

Pitch  cosine  component 
Pitch  Sine  Component 
Roll  Cosine  Component 
Roll  Sine  Component 
Antenna  Rate  Azimuth 
Antenna  Rate  Elevation 
Range  Fine 
Range  Coarse 
Spare 

Range  Rate  Fine 

True  Airspeed 

Vertical  Acceleration 

Log  of  Pressure  Coarse 

Spare 

Spare 

Log  of  Pressure  Fine 
Cosine  Antenna  Azimuth 
Angle 

Sine  Antenna  Azimuth  Angie 
Cosine  Antenna  Elevation 
Angie 

Sine  Antenna  Elevation 


INS  Synchro  Receive: 

INS  Synchro  Receive: 

INS  synchro  Receive: 

INS  synchro  Receive: 

Radar 

Radar 

Radar 

Radar 

Radar 

CADC 

INS 

CADC 


CADC 

Radar  Antenna 

Radar  Antenna 
Radar  Antenna 

Radar  Antenna 


TABLE  III-2  (sheet 


or  2  ) 


Measurement  Storage  in  Computer  '•lemorv 

DiMA 

LOCATION 

DESCRIPTION 

INPUT  SOURCE 

(OCTAL) 

040 

Range  Rate  Coarse 

Radar 

041 

Aspect  switch 

Radar  Set  Control 

Unit 

042 

Spare 

043 

spare 

044 

Spare 

045 

Spare 

046 

spare 

047 

Spare 

050 

spare 

051 

spare 

052 

Spare 

053 

spare 

054 

Spare 

055 

spare 

056 

Yaw  Cosine  Component 

INS  synchro  Receiver 

057 

Yaw  Sine  Component 

INS  Synchro  Receiver 

060 

spare 

061 

Spare 

062 

Spare 

063 

Spare 

064 

Spare 

065 

Spare 

066 

Spare 

06  7 

Spare 

070 

spare 

071 

spare 

072 

Spare 

o-o 

Spare 

074 

A.C.  Reference  Voltage 

Computer  Power 

Supply 

075 

Angle  of  Attack 

CADC 

076 

Spare 

0  77 

Spare 

»■«  w  ^  •r^ 


These  measurements  form  a  fixed  set  of  information  available 
to  a  Kalman  filter.  In  the  next  section,  measurements 
selected  from  this  set  are  mathematically  modelled  in  terms 
of  the  state  variables  defined  by  the  dynamics  models  of 
Chapter  II. 


I 


Measurement  Model  Equations 

In  the  previous  section,  t‘ e  entire  list  of  available 
measurements  produced  by  the  existing  F-4E/G  sensor  group 
was  examined.  From  this  group,  six  radar  measurements  are 
now  selected  as  Kalman  filter  inputs.  These  six  are  range 
(R),  range  rate  (R),  sine  of  the  antenna  azimuth  gimbal 
angle  (sin  X  ),  sine  of  the  antenna  elevation  gimbal  angle 
(sin  ) ,  azimuth  component  of  antenna  spatial  rate 
and  elevation  component  of  antenna  spatial  rate  (aj^).  The 
fine  scaling  for  R  and  R  is  used  since  this  is  the  normal 
scale  for  airborne  targets  inside  50  miles  (21:2-38).  In 
the  case  of  antenna  gimbal  angles,  the  cosine  information  is 
discarded  since  the  sine  completely  specifies  any  angle  in 
the  allowable  range  of  ±  60°.  This  section  develops  equa¬ 
tions  to  model  each  of  the  above  measurements. 

In  terms  of  the  position  states  defined  in  Fig.  2.2, 
the  range  from  the  radar  to  the  target  is  given  by 

R  =  (x^(t)  +  x^(t)  +  x2(t))^  (3-1) 


where  N,  E,  and  D  can  be  replaced  by  1,  3,  and  5  for  the  CV 
model  or  I,  4,  and  7  for  the  GMA  model.  Taking  the  time 


'A' 


f.'  j  ,  r: I"  I  m  P.'  pj  <■.1^. 


deriviative  yields  the  expression  f 

for  compactness,  R  is  expressed  in  terms  of 


time  arguments 


CV  states  as 

•  •  • 

*1*1  *3*3  *5*5 


R  = 


TJ 


,2  2  z, 

(x^  +  *3  +  X5) 


(3-2) 


Substituting  (2-3),  (2-5),  (2-7),  and  (3-1)  yields 


R  = 


*1^*2  '^FN^  *3(^4  “  '^FE^  *5^*6 


(3-3) 


For  mathematical  compatibility,  the  ownship  velocity  terms 
in  (3-3)  must  be  expressed  in  the  same  coordinate  system  as 
the  state  vector.  Thus,  \^p  must  be  expressed  in  the  geo¬ 
graphic  (N,E,D)  inertial  reference.  In  this  study,  super¬ 
script  G  implies  quantities  that  are  expressed  in  this 
specific  inertial  frame.  Initially,  the  F-4E/G  velocity 
vector  is  available  to  the  OFF  in  radar  reference  coordi¬ 
nates  (superscript  R),  as  shown  in  Appendix  C, 
transformation  1: 


TAS  cos(ai) 
0 

TAS  sin(Qt) 


(3-4) 


Obtaining  the  correct  velocity  components  requires  a  coord¬ 
inate  transformation  from  radar  reference  (i  ,i  ,k  )  coord- 

0  0  o 

inates  to  geographic  (N,E,D)  coordinates  given  by 

Ip  =  [T^]  (3-5) 
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where  [ Tp ]  is  the  required  transformation  matrix  (see 
K 

Appendix  C,  transformation  2). 

There  is  some  ambiguity  in  the  available  F-4E/G  docu¬ 
mentation  concerning  the  angles  and  X^ .  This  study 
assumes  that  they  are  the  actual  rauar  gimbal  angles  illus¬ 
trated  in  Fig.  3.3.  Azimuth  gimbal  angle,  X  ,  is  defined  as 
the  angle  between  the  radar  reference  i^-axis  and  the  projec¬ 
tion  of  the  line-of-sight  (los)  vector  on  the  i^'j^  plane. 
Elevation  gimbal  angle,  X^,  is  the  angle  between  the  los  and 
its  projection  on  the  Xq'Jq  plane.  Antenna  pointing  error 
is  ignored,  so  the  antenna  boresight  line  coincides  with  the 
los.  These  definitions  are  consistent  with  References  10 
and  20.  However,  Reference  14  describes  azimuth  and  eleva¬ 
tion  space  angles,  A^  and  E^ ,  used  in  the  fire  control 
system.  These  angles  are  defined  in  space  (x,y,z)  coordi¬ 
nates  as  shown  in  Fig.  3.4.  Aircraft  pitch  (0)  and  roll  ((J> ) 
relate  space  coordinates  to  aircraft  body  (l,m,n)  coordi¬ 
nates.  A  two  degree  offset  angle  between  the  body  1-axis 

and  the  radar  reference  i  -axis  is  described  in  Reference 

o 

10,  but  is  left  ambiguous  in  the  discussion  of  A^  and  E^  in 
Reference  14.  This  study  includes  the  two  degree  angle  in 
coordinate  transformations,  as  shown  in  Appendix  C. 

The  A^  and  E^  signals  are  used  exclusively  to  generate 
space  stabilized  cockpit  displays.  Pitch  and  roll  signals 
from  the  INS  are  sent  to  the  gyro-stabilized  platform 
(LRU  4)  where  they  are  used  to  convert  X  and  X  to  A  and 

3  S  S 

E  .  A  and  E  are  then  sent  to  the  control  indicator 
s  s  s 
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Target 


los  vector  can  be  used  in  right  triangle  relationships  to 


express  the  sine  of  either  gimbal  angle.  Next,  the  i^, 
and  components  are  expressed  in  terms  of  CV  or  GMA  state 

p 

variables  by  applying  the  inverse  of  the  [T^]  coordinate 
transformation  used  earlier.  In  radar  reference  coordi¬ 
nates,  the  trigometric  relation  for  the  sine  of  the  antenna 
azimuth  angle  is  given  by 


sin  A 


X  . 

- r- 

9  2 


(x“  +  x"  ) 

10  jo^ 


(3-6) 


The  required  transformation  is  illustrated  in  Fig.  3.5  and 
is  expressed  as 


x . 

x„ 

10 

N 

X  . 

JO 

=  (T«l 

’‘e 

1 

0 

X 

_ 1 

1^ 

X 

0 

1 

(3-7) 


Replacing  x^^  and  x^^  in  (3-6)  with  N,  E,  and  D  terms 
obtained  in  (3-7)  yields 


sin  A  = 
a 


^21’'n 


+  T22X£  +  ^23X5 


+  (T2  jXj^+T22X£+T23Xd) 


where  is  the  m,n-th  element  of  the  transformation  matrix 

described  in  Appendix  C,  transformation  3.  Substituting  the 
indices  1,  3,  5  or  I,  4,  7  for  N,  E,  D  gives  the  respective 
CV  or  GMA  expressions. 
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regardless  of  coordinate  frame.  The  previously  defined 
range  equation,  (3-1),  defines  this  quantity. 

The  remaining  radar  measurements  are  the  outputs  of  two 
antenna-mounted  rate  gyros.  The  gyros  measure  two  orthog¬ 
onal  components  of  the  antenna's  angular  rate  in  inertial 
space,  or  spatial  rate.  Expressed  in  antenna  (i,j,k)  coord¬ 
inates,  the  total  antenna  rate  is  represented  by  the  vector, 

A 

,  where  subscript  A  denotes  antenna  spatial  rate  and  super¬ 
script  A  indicates  the  rates  are  expressed  in  antenna  coord¬ 
inates.  One  rate  gyro  is  mounted  so  that  it  measures  the 
component  of  about  the  k-axis,  (Uj^,  and  the  other  gyro 
measures  the  component  of  about  the  j-axis,  . 

Fig.  3.6  illustrates  these  two  measurements.  Note  that  the 
third  component,  ,  is  meaningless  for  a  point  target  lying 
on  Che  i-axis.  The  target  will  always  lie  on  the  antenna 

i-axis  under  the  assumption  of  negligible  pointing  error, 

A 

and  thus,  oj  ,  and  uj  completely  characterize  w .  .  Since  the 
j  K  A 

above  assumption  implies  that  the  antenna  i-axis  and  the 

target  los  are  coincident,  the  measured  antenna  rate  also 

approximates  the  instantaneous  los  rate.  This  approximation 

can  be  used  to  model  the  antenna  rate  measurements. 

Using  u)  ,  and  u)  as  approximate  los  rate  components,  the 
J 

target's  relative  velocity  vector  can  be  expressed  in 
antenna  (i,j,k)  coordinates  by 


Note  that  an  equation  for  range  rate  also  results. 


R  '^FE^  ^13^ ’'6~'^FD^  (3-15) 


This  equation  could  conceivably  replace  the  range  rate  equa¬ 
tion  expressed  by  (3-3),  but  is  not  used  in  this  study.  The 
previously  described  range  rate  equation  is  superior  for 
practical  implementation,  since  [T]  matrix  elements,  a 
potential  source  of  error  in  finite  wordlength  computers  are 
not  required.  A  comparison  of  (3-3)  and  (3-15)  also  reveals 
an  alternate  method  of  calculating  the  first  row  of  the  [T] 
matrix,  by  equating  coefficients  of  common  terms. 

The  complete  measurement  models  for  both  CV  and  GMA 
filters  are  now  summarized.  The  selected  measurements  have 
been  deterministically  .xpressed  as  functions  of  both  the 
state  vector  and  independent  functions  of  time.  The  six 
measurements  are  now  combined  into  a  single  vector  of  func¬ 
tions  denoted  by  ti[2^(t),t]  and  defined  as 


This  definition  is  not  adequate  because  it  fails  to  describe 
the  effect  of  the  random  environmental  disturbances  which 


corrupt  all  physical  measurements.  To  account  for  these 
uncertainties,  additive  noise  is  included  in  the  measurement 
model.  The  vector  of  sampled-data  measurements  can  then  be 
expressed  as 


_z(t^)  =  il[i(t^),t^]  +  (3-17) 

where  t  ^ )  is  a  vector  of  zero-mean,  white  Gaussian  noise 
included  to  account  for  the  random  disturbances  which 
corrupt  the  measurements.  Substituting  the  expressions 
previously  derived  into  (3-17)  and  using  CV  state  variables 
yields  the  following  measurement  equations,  where  time  argu¬ 
ments  are  implicit: 


2  2 
+  X^  +  Xj 


i 

3 


(3-18) 


+  v , 


(3-19) 


2  .  . 2  .  ,2.^ 


(x^  +  x^  +  X.) 


^21^1  ^22''‘3  ^2  3''‘5 


( Ti  ^ X  i+T  1 2X3+T1 3X^  )  +  (Tt  jX^+T22X2  +  T-,2’‘5) 


n~  +  '^3 


(3-20) 


(T^iX^  +  T32’‘3  ^33^5) 

_  j 


+ 


(3-21  ) 
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he  tinal  consideration  needed  to  complete  the  Kalman 
r.lter  measurement  model  is  a  mathematical  characterization 
it  measurement  noise.  As  mentioned  earlier,  all  real  meas- 
'rements  of  physical  quantities  include  noise,  which  is  the 
combined  result  of  several  different  physical  effects.  Al¬ 
though  each  separate  effect  may  not  individually  possess  a 
(idussian  distribution,  the  combined  effect  of  several  small 
scale  effects  can  be  modelled  in  terms  of  Gaussian  random 
variables  by  virtue  of  the  central  limit  theorem.  It  is  the 


embodies  the  above  properties  in  the  statistics 


E{  v(t.)}  =  0  (3-24) 

E{  v(t.  )  J(tj)}  =  [R]  6.J  (3-25) 

where  E{  *}  is  the  expectation  operation  and  [R]  is  an  m  x  m 
matrix.  The  diagonal  terms  of  [R]  correspond  to  the  noise 
variances  of  the  six  measurements  in  the  measurement  vector. 
In  the  case  of  independent  measurements,  the  off-diagonal 
elements  of  [R]  are  all  zero.  However,  any  measurements 
which  exhibit  correlated  noise  properties  as  a  result  of 
physical  interrelationship  produce  non-zero  cross  terms  in 
[R]  matrix  row/column  intersections  common  to  both  measure¬ 
ments.  The  range  and  range  rate  measurements  produced  by 
the  F-4E/G  fire  control  system  are  examples  of  such  cross 
correlated  measurements.  They  are  both  processed  by  a  sin¬ 
gle  pair  of  cascade  integrators  in  the  electrical  synchro¬ 
nizer  (LRU  17),  and  share  the  thermal  noise  level  present  in 
that  circuit.  The  noise  correlation  is  not  total  because 
noise  sources  other  than  circuit  thermal  noise  independently 
corrupt  these  measurements. 

Actual  values  of  [R]  matrix  terms  can  be  determined  by 
theoretical  analysis  of  Che  sensor  devices  or  by  statistical 
reduction  of  experimental  data.  Both  techniques  are  used  in 
this  study.  Engineers  at  Hill,  AFB,  Utah  provided  recom¬ 
mended  variances  for  range,  range  rate,  gimbal  angle,  and 
angle  rate  measurements  based  on  actual  flight  tests  (24). 


This  information  is  augmented  by  data  from  a  computer  ex¬ 
periment  used  to  determine  variances  for  sin  and  sin 
Typical  values  of  conical  scan  AI  radar  accuracies  listed  in 
standard  textbooks  and  F-4E/G  maintenance  manuals  are  also 
compared  during  filter  tuning  (14,  21,  25).  Table  III-3 
summarizes  the  actual  variance  values  used  for  both  the  CV 
and  GMA  measurement  models. 

TABLE  III-3 


Measurement  Noise  Variances  for  F-4E/G  LRI  Filter 


Measurement 

l-Sigma 

Variance 

Range 

67  feet 

4489 

Range  Rate 

16  f ps 

256 

Sin  X 

a 

1.05  X  10"^ 

1.10  X  10"^ 

Sin  X 

e 

1.05  X  10“^ 

1.10  X  10"^ 

0)  . 

1 

3.49  X  10  ^  rad/ s 

1.22  X  10"^ 

3.49  X  10  ^  rad/s 

1.22  X  10~^ 

An  additional  consideration  in  determining  appropriate 
[R]  matrix  values  for  this  study  is  the  constrained  word- 
length  of  OFF  measurements.  Although  the  LRU-1  computer 
itself  has  16-bit  wordlength,  the  A/D  converter  which  pro¬ 
cesses  the  measurements  prior  to  storage  in  the  computer 
allows  only  10-bit  wordlength.  This  restriction  can  affect 
the  variance  of  measurement  error  under  certain  conditions. 


•.^1"  ^ ^  ■ '-  ■'.'  ■  ■  '1 
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For  example,  at  maximum  LRI  ranges,  the  A/D  converter's 
wordlength  restriction  causes  the  least  significant  bit  of 
the  range  measurement  to  have  a  value  of  256  feet.  This  is 
effectively  a  3-sigma  value  for  the  range  measurement  and 
corresponds  to  a  variance  of  7282,  or  about  twice  the  range 
variance  listed  in  Table  III-3. 

The  above  quantification  of  measurement  noise  completes 
the  development  of  models  needed  to  form  the  basis  of  a 
Kalman  filter.  The  implementation  of  these  models  in  an 
extended  Kalman  filter  format  is  the  subject  of  the  next 
chapter . 


•  •  -  •  I 
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IV .  Filter  Implementation 


Extended  Kalman  Filter 

Model  nonlinearities  are  an  intrinsic  characteristic  of 
the  air-to-air  tracking  problem.  As  discussed  in  Chapter  II 
the  state  variables  of  the  CV,  GMA,  and  other  dynamics 
models  can  be  defined  as  total  or  relative  states  and  refer¬ 
enced  to  any  of  several  inertial  or  line-of-sight  coordinate 
frames.  Regardless  of  the  state  variable  definitions  and 
reference  frame  used,  each  formulation  leads  to  nonlineari¬ 
ties  in  either  the  dynamics  model,  the  measurement  model,  or 
both  (15,16).  It  is  these  inherent  nonlinearities  of  the 
tracking  problem  which  preclude  exclusive  use  of  conven¬ 
tional  linear  Kalman  filter  techniques  in  this  study.  Pos¬ 
sible  substitutes  include  linearized  Kalman  filters,  ex¬ 
tended  Kalman  filters,  second-order  filters,  or  other  forms 
of  state  estimators  such  as  maximum  likelihood  estimators 
(2:24,212).  The  extended  Kalman  filter  (EKF)  is  selected  as 
the  most  suitable  implementation  based  on  the  three  reasons 
explained  in  the  introductory  chapter:  proven  performance, 
computational  t r ac t a b i 1 i t y ,  and  compatibility  with  available 
evaluation  tools.  The  enhanced  tractablility  of  the  EKF 
results  from  linearizing  the  models  about  nominal  points  on 
the  estimated  state  trajectory,  although  some  portions  of 
the  algorithm  still  require  the  explicit  evaluation  of  non- 
linearities.  The  EKF  maintains  accuracy  by  repeatedly 


relinearizing  the  filter  about  each  new  state  estimate,  in 
effect  providing  a  new  nominal  trajectory  for  each  succes¬ 
sive  filter  operation. 

This  chapter  explains  EKF  formulations  for  the  models 
developed  in  preceding  chapters.  Filter  propagation  between 
measurement  times  is  developed  first.  Both  extended  and 
conventional  formulations  are  discussed  in  the  filter  propa¬ 
gation  section  because  the  preservation  of  linear  dynamics 
in  the  CV  and  GMA  filters  allows  the  exploitation  of  linear 
techniques  to  simplify  propagation.  The  next  section  dis¬ 
cusses  the  EKF  measurement  update  equations  required  by 
nonlinearities  in  the  measurement  models  of  both  filters. 

The  chapter's  final  section  describes  an  EKF  algorithm  for 
the  F-4E/G  OFF  and  the  evaluation  of  several  key  parameters. 
The  actual  EKF  equations  and  their  equivalent  discrete  ver¬ 
sions  are  listed  without  proof.  A  detailed  derivation  is 
available  in  Reference  2,  Chapter  9. 

Filter  Propagation 

In  the  current  context,  the  term  propagation  refers  to 
the  manner  in  which  the  states  of  a  given  process  change,  or 
evolve,  in  time.  Applied  to  the  sampled  data  Kalman  esti¬ 
mator,  propagation  provides  a  description  of  the  time  evolu¬ 
tion  of  the  states  during  the  periods  of  time  between  meas¬ 
urements.  A  Kalman  filter  propagates  both  the  vector  of 
state  estimates,  i(t^),  and  its  associated  covariance 


matrix ,  P( t . ) . 


The  Kalman  filter's  internal  dynamics  model  governs 
filter  propagation.  For  the  extended  Kalman  filter  formula 
tion,  system  dynamics  are  modelled  in  the  general  form 

2c(t)  =  f^[x^(t)  ,^(t)  ,t  ]  +  G(t)w(t)  (4-1 

where 

jc  ( t )  is  the  state  vector, 

f^[  Jt  ( t )  ,^(  t )  ,  t  ]  is  a  vector  of  functions  describing 
state  dynamics  (not  necessarily  linear), 

^(  t )  is  a  vector  of  deterministic  inputs, 

G(t)  is  a  coefficient  matrix  defining  which  states  are 
driven  by  linearly  additive  noise  terms, 

and 

^(t)  is  a  vector  of  zero-mean  white  Gaussian  noise  of 

strength  Q(t)  as  defined  by  (2-9). 

Defining  an  initial  condition  for  (4-1)  as  a  random  vector 

x(t  )  with  mean  x  and  covariance  P  completes  the  time 
~  o  — o  o 

domain  description  of  the  system  of  interest. 

Given  discrete  time  measurements  up  to  and  including 
time  t^,  the  propagated  state  estimate  at  any  subsequent 
time  t  prior  to  the  next  sample  time  is  denoted  £(t/t^)  and 
is  obtained  by  integrating  the  equation 

£(t/t^)  =  f_[£(t/t^)  ,u^(t),t]  (4-2 

forward  to  time  t  from  the  initial  condition 

x(t./t.)  =  x(t^)  (4-3 

where  i(tt)  describes  the  result  of  the  previous 


measurement  update  at  time  t^.  Superscript  +  indicates 
estimates  immediately  after  the  incorporation  of  measure¬ 
ments.  If  the  recursion  begins  with  a  propagation,  the 
initial  condition  for  the  first  propagation  is 

x(t/t)=x  ( 4-4) 

—  00—0 

The  propagated  error  covariance  matrix  of  the  estimated 
state  vector  is  calculated  similarly.  Integrating  the  equation 

P(t/t.)  =  F[t;x(t/t  .  )  ]P(t/t.  )  +  P(t/t  .  )F'^[t;x(t/t  .  )  ] 

+  G(t)Q(t)G'’’(t)  (4-5) 

forward  to  time  t  from  the  initial  condition 

P(t./t.)  =  P(t+)  (4-6) 

yields  the  newly  propagated  state  covariance  matrix 
P(t“^^),  where  superscripts  -  and  +  respectively  denote 
quantities  before  and  after  measurement  update.  The  ini¬ 
tial  condition,  defined  above  by  (4-6),  is  obtained  from 
the  filter's  previous  measurement  update  or  from  the  val¬ 
ues  in  P  in  the  case  of  the  initial  recursion.  The 
o 

matrix  F [  t  ;_x ( t / 1 ^  )  ]  in  (4-5)  represents  the  Jacobian  of 
the  newly  relinearized  dynamics  model  and  is  defined  by 

t ) ,  t  ] 

F[t;i(t/t  )]  =  -  (4-7) 

9)^  3^=:£(t/t^) 

where  i(t/t^)  is  the  current  state  estimate  obtained  by 
integrating  (4-1). 


For  computer  implementations,  a  discrete-time  form  of 


filter  propagation  is  required.  In  this  form,  a  newly 
propagated  state  estimate  for  time  is  obtained  by 
integrating  (4-1)  to  yield: 


x(t.  ,)  =  x(t.  ,/t.) 
—  1  +  1^  —  1+1  1 


(4-8) 


This  result  is  used  to  evaluate  the  current  F  matrix,  ex¬ 
pressed  in  (4-7),  which  is  needed  in  the  covariance  propaga¬ 
tion  equation,  defined  by  (4-5).  The  relinearized  covar¬ 
iance  propagation  equation  is  integrated  next  to  yield  the 
propagated  state  covariance  matrix  at  time 


This  completes  the  EKF  propagation  from  t^  to 

If  the  function  f_[£(  t  )  ,^(  t ) ,  t  ]  is  nonlinear,  as 
implied  by  a  nonlinear  dynamics  model,  solving  the  propaga¬ 
tion  equations,  (4-1)  and  (4-5),  requires  some  form  of 
numerical  integration  in  the  computer  implementation.  The 
computational  load  of  this  requirement  is  severe.  However, 
the  propagation  equations  can  be  simplified  if  system  dy¬ 
namics  are  modelled  using  linear,  time-invariant  differen¬ 
tial  equations.  If  the  dynamics  model  is  linear,  (4-1)  can 
be  expressed  in  matrix  form  as 


x(  t ) 


F(t)£(t)  +  B(t)u(t)  +  G(t)w(t) 


(4-10) 


where  F(t)  is  the  plant  matrix  and  B(t)  is  called  the  input 
matrix.  The  solution  to  (4-10)  is  given  by 
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where  ^  state  transition  matrix  which  propa¬ 
gates  the  solution  from  time  t^  to  time  snd 

the  initial  condition  for  (4-10).  The  second  integral  term 
in  (4-11)  is  a  stochastic  integral  representing  discrete¬ 
time  white  Gaussian  noise,  derived  from  a  Brownian  motion 
process,  §_( t ) .  Reference  1,  Chapter  4  explains  the  details 
of  this  derivation.  If  the  system  is  time-invariant,  F(t) 
is  a  constant  matrix  and  ( t .  ..t.)  is  a  function  of  the 
time  difference,  the  propagation  interval 

is  constant,  as  is  the  case  when  the  time  between  measure¬ 
ments  is  equal,  the  state  transition  matrix  is  a  constant 
matrix  which  can  be  precomputed  for  all  time.  A  Kalman 
filter  operating  under  the  above  conditions  propagates  state 
and  covariance  estimates  from  time  tt  to  time  t.  ,  by  the 
equations 


_x(t^^j)  =  <{)  (At  )2L(  t^  )  +  B^(t^)u(t^) 
P(t~^P  =  <I>(At)P(t^)<t'^(At)  +  Qd(t.) 

where 

=  t ,  .  ,  -  t ,  , 


(4-12) 

(4-13) 


At 


(4-14) 


(4-15) 


j 


<l>(t 

t  . 

1 


i  +  1 


.T  )B(T  )dT  , 


’^i+l 


Qd(t. 


.)  =  y  <I>(t^^^,T)G(T)Q(T)G’^(T)<l>'^(t.^^,T)dT 


(4-16) 


and  ii(t^)  is  an  input  vector  approximated  as  constant  over 
each  sample  period.  The  constant  state  transition  matrix 
<l>(At)  can  be  evaluated  in  closed  form  using  Laplace  trans¬ 
form  theory: 


<I>(At) 


L  ^[sI-F]  ^ 


At 


(4-17) 


A  filter  propagation  scheme  which  employs  (4-12)  and  (4-13) 
can  be  implemented  computationally  by  straightforward  matrix 
addition  and  multiplication.  Herein  lies  the  great  advan¬ 
tage  of  preserving  linear,  time-invariant  dynamics  models 
and  incorporating  measurements  taken  at  equal  intervals. 
Propagation  under  these  conditions  avoids  the  EKF  require¬ 
ments  for  the  computationally  expensive  numerical  integra¬ 
tion  and  recursive  evaluation  of  F [ t ;^( t / t ^ ) ] . 

As  discussed  in  Chapter  II,  linear,  time-invariant 
dynamics  are  intentionally  preserved  in  the  CV  and  GMA 
filter  designs.  Additionally,  the  APQ-120  takes  measure¬ 
ments  at  a  constant  rate,  every  0.01  second,  and  the  OFP  can 
be  programmed  to  incorporate  measurements  at  any  multiple  of 
this  period  (8).  Thus,  both  models  satisfy  the  prerequi¬ 
sites  for  propagation  by  means  of  (4-12)  and  (4-13). 


For  the  CV  model,  putting  the  dynamics  equations, 
(2-3)  through  (2-8),  into  the  form  of  (4-10)  results  in 


Calculating  the  state  transition  matrix  for  At  =  0.04 


seconds  yields 


1  0.04  0  0  0  0 

0  1  0  0  0  0 

<f(At)  =  0  0  1  0.04  0  0  (4-19) 

0  0  0  1  0  0 

00001  0.04 

_  0  0  0  0  0  1  _ 


as  shown  in  Appendix  B. 


n»  ■  "v  • 

K 


*  ^ 

*  •  • 

i:. 

» , 


Evaluating  (4-15)  and  (4-16)  for  the  0.04  measurement 
interval  yields  the  discrete-time  input  and  driving  noise 


matrices , 


"d  = 


-0.04  0 


0  0 

-0.04  0 

0  0 

0  -0.04 

0  0 


(4-20) 


^d  = 


0  0 

0  0.04q 

0  0  ' 

0  0 

0  0 

0  0 


0  0.04q  0 
0  0  0 


0  0 

0  0 

0  0 

0  0 

0  0 

0  0.04q, 


(4-21) 


where  q^^,  q^,,  and  q^  in  (4-21)  are  constant  noise  strengths 

associated  respectively  with  w(t),  w(t),  and  w(t),  the 

N  c  D 

elements  of  the  driving  noise  vector. 

The  determination  of  $(At),  B ,  ,  and  Q,  completes  the 

Q  u 

linear,  time-invariant  version  of  the  recursive  propagation 
expressions,  (4-12)  and  (4-13),  as  applied  to  the  CV  model. 
To  begin  the  recursion,  the  filter  requires  initial  condi¬ 
tions  X  and  P  .  For  this  study,  the  initial  state 
estimate  is 


X  =  0 

— o  — 


(4-22) 


Values  for  the  initial  covariance  matrix  P  are  first  se- 

o 

lected  as  reasonable  estimates  of  the  filter's  initial  state 
uncertainty.  However,  these  values  can  be  adjusted  during 
the  filter  tuning  process  to  improve  the  initial  convergence 
performance  of  the  filters. 


The  propagation  equations  for  the  GMA  model  are  found 


by  a  similar  process.  Modelling  the  dynamics  equations, 
(2-11)  through  (2-19),  in  the  form  of  (4-10)  yields 


(  t  ) 

X  j  ( t ) 

X2(  t ) 

X2(  t ) 

^^(t) 

X3(  t ) 

I_  — . 

l_  — , 

^^(t) 

X4(t) 

w^(t) 

X2(  t  ) 

=  [F] 

Xg  ( t ) 

+  [B] 

+  [G] 

WE(t) 

X^(t) 

x^(t) 

WD(t) 

x^(t) 

XyCt) 

igCt) 

X^it) 

X9(  t ) 

X9(t) 

(4-23) 

where 
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The  input  and  noise  coefficient  matrices  for  the  GMA  model 
are 


-1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

-1 

0 

0 

0 

0 

0 

0 

0 

and 

G  = 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

-1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

In  Appendix  B,  the  state  transition  matrix  is  calculated 


<t>(At) 


where 


and  At  = 

0.04  seconds: 

3x3 

1  °3x3 

_°3x3 

_°3x3 

1  * 

3x3 

1 

_‘^3x3 

(4-25) 

°3x3 

^  °3x3 

1 

*3,3^ 

~1  0.04 

0.0007791“ 

0  1 

0.0384418 

0  0 

0.9231163 

calculated  using  (4 

-15)  and 

(4-16) 

.  Initial 

3x3 


d  "d 

conditions  x  and  P  are  chosen  in  the  same  manner  as  in 


— o  o 

the  CV  model,  but  with  additional  terms  for  the  acceleration 
states.  Using  't(At),  and  just  described, 

the  GMA  filter  can  be  propagated  by  the  linear  Kalman  filter 
equations,  (4-12)  and  (4-13). 

This  completes  the  description  of  time  propagation 


in  the  CV  and  GMA  filters 


In  these  models,  linear 


time-invariant  dynamics  lead  to  simplified  calculation  of 
state  and  covariance  estimates,  providing  the  estimates  are 
propagated  over  equal  time  intervals  between  measurements. 
The  next  section  describes  the  EKF  equations  for  incorpo¬ 
rating  measurements  into  the  state  estimate  at  each  measure 

ment  time  t . . 

1 

Measurement  Update  for  Extended  Kalman  Filters 

At  each  measurement  time  t.,  the  EKF  receives  informa- 

1 

tion  from  sensors  in  the  form  of  a  discrete  time  measurement 
vector  Chapter  III  developed  a  mathematical  model  of 

the  measurements  available  to  the  F-4E/G  fire  control  com¬ 
puter.  This  model  expresses  as  the  sum  of  nonlinear 

functions  of  state  variables  and  zero-mean  white  Gaussian 
noises : 


£(t^)  =  h[x(t^),t^]  +  v(t^)  (4-26) 

which  is  in  the  form  of  (3-17).  This  model  governs  the 
EKF  update.  For  each  update  cycle,  the  filter  algorithm 
evaluates  the  measurement  function  vector  h.[*,*]  using  the 

A  _ 

current  propagated  state  estimate  x(ti).  This  represents 

the  filter's  best  prediction  of  what  each  measurement  should 

be  at  time  t . : 

1 

£(tp  =  1l[  x.(  t”)  ,  t^  ]  (4-27) 

When  the  actual  noise-corrupted  measurements  from  the  en¬ 


vironment,  z.,  become  available,  the  EKF  calculates  the 


difference  between  the  actual  and  predicted  values.  The 
difference,  or  residual,  is  multiplied  by  an  optimal  weight¬ 
ing  matrix  K(t^)  to  produce  a  correction  term  for  the  state 
estimate.  In  this  manner,  the  EKF  incorporates  measurement 

A  _ 

information  into  the  current  propagated  state  estimate  x^t  | 
to  produce  an  updated  state  estimate  for  the  current  time, 
^(t^).  In  equation  form  this  is 

i(tt)  =  i(t^)  +  K(t^){z^^  -  h[^(tj),t^]}  (4-28) 

where  is  the  vector  of  actual  measurements.  The  optimal 
weighting  matrix  K(t^)  is  found  by 

K(t.)  =  P(t^)H^(t.  )  {H(t.  )P(tpH'^(t.  )  +  (4-29) 

where  R(t^)  is  the  measurement  noise  matrix  described  in 
Chapter  III  and  H(t^)  is  a  partial  deriviative  matrix 
defined  by 

9  h.  [21,  t .  ] 

H[t  ;B(t~)]  =  - ^  ^  _  (4-30) 

9_X  3C  =  ^(  t^  ) 

The  covariance  is  updated  by 

=  P(t“)  -  K(t.  )H(t.  )P(t]^)  (4-31) 

where  K(t^)  and  H(t^)  are  defined  respectively  by  (4-29)  and 
(4-30) . 

The  measurement  function  vector  h_[*,*],  as  defined  by 
(3-16),  is  physically  the  same  for  both  the  CV  and  GMA 
models.  As  illustrated  below,  this  results  in  mathematical 


similarities  in  the  H  matrix  of  both  models.  The  only 
differences  are  attributable  to  the  different  dimensionality 
and  the  state  variable  numbering  conventions  of  the  two 


models.  Calculating  the  partial  deriviative  matrix  indi¬ 
cated  by  (4-30)  in  terms  of  the  CV  state  vector  yields 


«1 

0 

«2 

0 

H3 

0 

»4 

«5 
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«7 
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»9 

H  = 

^10 

0 

«11 
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«15 
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^18 

Hi9 
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^23 

»24 

H25 

»26 

H27 

where  individual  elements  are  defined  by  the  following 
equations.  Zero  valued  elements  are  not  listed. 


■1 


,  2  2  2,4 

(Xi  +  X3  +  X5) 


T2 2 2,4 

(x^  +  X3  +  x^)^ 


T2  2  2T4 

(X^  +  X3  +  X3) 


(A-33) 


(4-34) 


(4-35) 


(x; 


)  - 


XjX3(x^- 


) 

FE'’ 


(x: 


+  X. 
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V. 


V 


(4-36) 


M33  I  - 


(4-46) 


-TSsiCx^  +  X3)  +  1332^2^3  +  ^^33^1^5 
ri  2  2.3/2 

(Xj  +  X3  +  X3) 


(4-47) 


T332(Xi  +  X3)  +  1333X^X3  +  T333X3X3 

,  2  2  2,3/2 

(Xi  +  X3  +  X3) 


(4-48) 


-  T333(x^  +  X3)  +  T33j^Xj^X3  +  T332X3X3 

=  2  7-  2,3/2 

(X)  +  X3  +  X3) 


(4-49) 


Tiie  remaining  H  terms  require  [T4],  Transformation  4  in 
Appendix  C. 


Xj^[T43^(x2  Vpj^)  +  T432(x^  V^^,)  +  T433(x^  ^ 


,  1  ^  2  ^  2,3/2 

(X^  +  X3  +  X5) 


(4-50) 


,  2  2  2,^ 

(x^  +  X3  +  x^)^ 


(4-51) 


^3^T43i(x2  Vp^)  +  T^32^^4  '^fE^  ^^33^^6 


,  1  2  ^  2,3/2 

(X^  +  X3  +  X3) 


(4-52) 


,2  2  2,^ 

(X^  +  X3  +  X3) 


(4-53) 


(4-62) 


Algorithm  for  OFF  Implementation 

This  section  provides  a  step-by-step  procedure  which 
can  be  used  as  a  framework  for  programming  the  EKF  implemen¬ 
tations  discussed  in  this  chapter.  The  three  basic  steps 
are  filter  initialization,  propagation,  and  update.  These 
steps  are  expanded  below. 

(1)  The  following  vectors  and  matrices  are  initialized 
at  program  start. 

(a)  Noise  matrices.  Initial  driving  noise  strength 
Q(t),  is  discussed  in  Chapter  II.  The  equiv- 
ilent  discrete  noise  matrix,  ,  needed  for 
this  algorithm  is  given  by  (4-16),  with  an  ex¬ 
ample  calculation  in  Appendix  B.  Chapter  III 
discusses  measurement  noise  strengths.  Final 
noise  strengths  are  obtained  during  the  tun¬ 
ing  process,  as  described  in  Chapters  V  and  VI. 


(b)  State  transition  matrix.  This  matrix  is 
calculated  by  applying  (4-17)  to  a  linear, 
time-invariant  dynamics  model.  It  is  con¬ 
stant  for  a  given  measurement  interval. 
Appendix  B  contains  sample  calculations. 

(c)  Input  matrix.  is  calculated  using  (4-15) 

An  example  is  in  Appendix  B. 

(d)  Initial  conditions.  The  initial  state  esti¬ 
mate  and  covariance  matrix  are  discussed  in 
the  filter  propagation  section.  For  ,  the 
values  are  often  adjusted  during  filter  tun- 


ing . 
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attitude 

angles  are  needed  to 

calculate  transformation  matrices  and  veloc¬ 
ity  aiding  terms  for  the  initial  propagation 

(2)  Propagation. 

(a)  Calculate  the  propagated  state  estimate, 
using  (4-12). 


(b)  Calculate  the  propagated  covariance  matrix, 
using  ( 4- 1 J  )  . 


( 3  )  Update . 

When  each  new  set  of  measurement  data  becomes 
available,  the  following  steps  are  accomplished  in 
the  order  given: 

(a)  Calculate  transformation  matrix  elements 
based  on  the  new  aircraft  attitude  angles. 

(b)  Calculate  current  ownship  velocity  components 
in  north,  east,  and  down  coordinates. 

(c)  Calculate  the  predicted  measurement  vector, 
(4-27),  using  the  result  from  step  (2a)  above. 

(d)  Calculate  the  partial  deriviative  matrix, 

(4-30),  evaluated  at  ^(t^),  the  result  from 
step  ( 2a ) . 

(e)  Calculate  the  Kalman  filter  gain,  (4-29). 

(f)  Calculate  the  updated  state  estimate,  (4-28). 

(g)  Calculate  the  updated  covariance  matrix,  (4-31). 

(4)  Begin  the  next  propagation  cycle  and  continue  the 
recursion . 

The  above  algorithm  is  designed  to  implement  an  operational 
extended  Kalman  filter  in  the  F-4E/G  fire  control  system. 

A  method  of  evaluating  this  filter's  performance  without 
actually  implementing  it  on  an  aircraft  is  described  in  the 
next  chapter. 


V.  Method  of  Evaluation 


Introduction 


An  important  advantage  of  Kalman  filtering  is  its  suit- 


ablity  for  computer  testing  and  evaluation.  Computer  simu¬ 


lation  allows  repetitive  evaluation  of  filters  without  the 


expense  of  constructing  multiple  prototypes.  For  this  rea¬ 


son,  the  filters  developed  in  this  study  are  evaluated  in  a 


computer  simulation.  This  chapter  discusses  the  structure 


of  the  simulation  software  and  the  methods  employed  to 


evaluate  filter  performance  using  this  simulation.  The 


first  section  introduces  a  dedicated  software  package  de¬ 


signed  specifically  for  the  evaluation  of  extended  Kalman 


filters.  A  second  section  describes  a  computer  program 


acquired  from  F-4E/G  software  engineers  at  Hill  AFB,  Utah 


and  modified  to  provide  a  standard  of  comparison  for  filter 


performance.  This  standard,  or  truth  model,  augments  the 


basic  Kalman  filter  evaluation  software  to  form  a  complete 


simulation  package  for  testing  the  CV  and  GMA  filters.  The 


final  section  discusses  the  general  procedures  followed  in 


using  the  simulation  package  for  successive  evaluations. 


improvements,  and  re-evaluations  of  filter  performance. 


SOFE 


A  number  of  software  packages  currently  exist  for  the 


purpose  of  evaluating  Kalman  filter  designs.  One  such  pro¬ 


gram  is  SOFE,  a  general  purpose  Monte  Carlo  simulation,  and 


,v' •. 


a  companion  data  postprocessor  and  plotting  program,  SOFEPL 
(13,28).  Monte  Carlo  analysis  is  a  method  of  calculating  a 
statistical  description  of  a  stochastic  process  and  requires 
repeated  runs  of  a  computer  simulation  to  generate  many 
samples  of  the  stochastic  process.  Each  sample  varies  under 
the  influence  of  a  random  number  generator  used  as  the 
source  of  realized  pseudo-random  effects  within  the  process. 
Mean  and  variance  statistics  are  calculated  for  each  time 
these  statistics  are  needed  throughout  the  simulation  by 
averaging  data  generated  over  the  entire  set  of  samples.  If 
the  number  of  samples  is  sufficiently  large,  the  process 
statistics  from  the  Monte  Carlo  study  closely  approximate 
the  statistics  of  the  real  world  process  (1:329).  In  using 
the  Monte  Carlo  technique,  the  analyst  must  address  the 
issue  of  how  many  samples  are  sufficient  to  characterize  the 
process  statistics  adequately.  Previous  studies  suggest  20 
runs  as  adequate  for  many  applications  (12:65,27:59).  In 
this  study,  5  run  studies  are  used  to  examine  general 
trends,  and  20  run  studies  are  used  to  fully  characterize 
process  statistics  in  the  simulation. 

In  the  present  case,  the  stochastic  process  to  be 
statistically  characterized  is  a  vector  error  process  con¬ 
sisting  of  the  errors  committed  by  a  reduced-order  Kalman 
filter  as  it  estimates  target  position,  velocity,  and 
acceleration  values  in  the  F-4E/G  OFP.  The  filter's  true 
errors  are  the  differences  between  the  filter-computed 
values  and  the  true  values  of  target  position,  velocity,  and 


acceleration  in  a  real  world  scenario.  In  actual  air  com¬ 
bat,  the  true  values  of  these  quantities  are  not  known.  In  a 
computer  simulation,  however,  the  true  values  of  target 
position,  velocity,  and  acceleration  can  be  approximated  by 
a  mathematical  model.  Additionally,  models  can  approximate 
true  measurements  which  emulate  the  response  of  an  APQ-120 
radar  operating  under  conditions  associated  with  specific 
F-4E/G  combat  maneuvers.  A  mathematical  model  which  closely 
approximates  the  true  values  of  variables  of  interest  within 
a  dynamic  system  is  called  a  truth  model.  A  truth  model 
must  be  the  most  accurate  and  complete  model  available 
because  the  validity  of  the  filter's  performance  evaluation, 
and  ultimately,  the  operational  success  of  the  filter,  de¬ 
pend  on  the  adequacy  of  the  truth  model  to  represent  reality 
(1:326).  Using  a  suitable  truth  model,  a  sample  of  the 
filter's  error  process  can  be  generated  by  subtracting  the 
outputs  of  the  truth  model  from  the  corresponding  outputs  of 
the  filter  at  each  discrete  time  of  a  simulation  run.  Addi¬ 
tional  samples  are  generated  in  the  same  way,  but  remain 
independent  of  other  samples  through  the  use  of  a  non- 
repetitive  random  number  generator.  Each  sample  associates 
a  vector  of  error  states  with  each  discrete  time  of  the 
simulation.  For  each  of  these  discrete  times,  summing  the 
corresponding  components  of  the  error  state  vectors  asso¬ 
ciated  with  that  time  and  dividing  by  the  number  of  samples 
yields  a  mean  error  state  vector  for  that  particular  time. 
Similarly,  the  error  variance  for  each  discrete  time  is  also 


calculated  directly  from  the  data  generated  by  the  samples. 
When  these  calculations  are  done  for  all  times  in  the  simu¬ 
lation,  the  completed  Monte  Carlo  study  characterizes  the 
time  history  of  the  mean  and  variance  of  the  error  process 
associated  with  each  component  of  the  filter  state  vector. 
The  following  discussion  illustrates  Monte  Carlo  analysis  in 
the  context  of  the  SOFE  program. 

Since  SOFE  is  a  general  purpose  Monte  Carlo  simulation, 
it  can  be  adapted  to  a  variety  of  Kalman  filter  applica¬ 
tions.  The  user  defines  a  specific  problem  by  supplying 
both  the  filter  structure  and  the  truth  model.  User-speci¬ 
fied  input  variables  determine  the  number  of  filter  states, 
the  number  of  measurements,  the  measurement  time  increment, 
and  other  parameters  which  govern  the  dimensions  of  the 
models,  program  control,  and  input/output  timing.  The 
actual  filter  dynamics  and  measurement  equations  are  imple¬ 
mented  in  user-written  subroutines  appended  to  the  basic 
SOFE  software.  The  basic  SOFE  code  itself  provides  the 
overall  Monte  Carlo  structure,  control  for  event  timing, 
and  computational  tools  for  propagating  and  updating  a 
Kalman  filter  state  vector  and  its  associated  covariance 
matrix  (13). 

The  SOFE  user  provides  a  truth  model  in  one  of  two 
ways.  It  can  be  implemented  in  user-written  subroutines,  a 
method  which  allows  SOFE  to  propagate  the  truth  states  over 
each  time  step  in  parallel  with  the  filter  states,  or  alter¬ 
nately,  the  truth  model  can  be  supplied  by  means  of  an 


external  state  trajectory.  The  latter  configuration  re¬ 
quires  a  separate  data  file  that  contains  a  time  history  of 
truth  state  values  and/or  true  measurements.  This  file  can 
consist  of  actual  data  recorded  during  operational  tests  of 
the  physical  system.  Frequently,  however,  the  physical 
system  is  not  yet  built  or,  as  is  the  case  in  this  study, 
not  readily  available.  In  this  case,  a  separate  computer 
program  serves  as  the  truth  model  and  generates  the  external 
trajectory  data  a  priori  for  storage  in  the  data  file.  SOFE 
accesses  and  interpolates  data  from  this  file  as  necessary 
during  the  course  of  each  Monte  Carlo  run.  The  external 
trajectory  truth  model  configuration  is  particularly  well 
suited  to  target  tracking  problems  and  is  the  choice  for 
this  study  (29).  The  actual  program  used  is  covered  later 
in  this  chapter. 

The  basic  SOFE  code  includes  subroutines  which  propa¬ 
gate  and  update  both  truth  and  filter  states.  SOFE  accom¬ 
plishes  filter  propagation  by  means  of  a  fifth-order  Kutta- 
Merson  numerical  integration  routine.  The  simulation  propa¬ 
gates  filter  estimates  by  this  method  regardless  of  the 
linearity  or  nonlinearity  of  the  dynamics  model.  The  time 
interval  between  measurements  is  subdivided  into  adjustable 
length  subintervals.  SOFE  integrates  the  propagation  equa¬ 
tions  for  both  the  state  vector  and  covariance  matrix  over 
each  subinterval  using  the  results  from  the  previous  sub¬ 
interval  integration  as  initial  conditions.  Also,  the  state 
estimate  obtained  by  integrating  over  the  previous  subinterval 
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is  used  to  evaluate  the  F  matrix,  defined  by  (4-7),  which  is 
needed  to  propagate  the  covariance  over  the  current  subin¬ 
terval.  SOFE  automatically  adjusts  the  size  of  each  subin¬ 
terval  to  maintain  integration  errors  within  user-defined 
tolerances  while  minimizing  computer  execution  time 
(30:747).  As  propagation  progresses,  SOFE  monitors  the 
current  time  and  compares  it  to  the  next  scheduled  measure¬ 
ment  time.  When  a  calculated  subinterval  reaches  or  exceeds 
the  next  measurement  time,  SOFE  truncates  that  subinterval 
at  the  measurement  time  and  terminates  the  propagation  cycle 
at  the  end  of  that  subinterval.  At  that  point,  SOFE  updates 
the  filter  state  estimates  and  covariance  by  means  of  a 
sequential  Carlson  square  root  algorithm,  an  update  method 
with  well  established  numerical  stability  (1:385,13:26). 

The  Carlson  method  allows  vector  measurements  to  be  incorpo¬ 
rated  sequentially  as  scalars.  This  capability  allows  SOFE 
users  to  omit  some  measurements  during  certain  Monte  Carlo 
studies  in  order  to  determine  the  effect  of  specific  meas¬ 
urements  on  filter  performance.  When  all  desired  measure¬ 
ment  updates  are  complete,  SOFE  resumes  the  propagation 
cycle.  Propagation  and  update  cycles  continue  to  alternate 
until  the  end  of  each  Monte  Carlo  run. 

As  SOFE  accomplishes  each  Monte  Carlo  run,  it  stores 
truth  states,  filter  states,  filter  variances,  residuals, 
and  residual  variances  for  each  incremental  filter  time  on 
an  external  interface  file  (28:5).  Following  the  completion 
of  all  runs,  the  postprocessing  and  plotting  program, 


SOFEPL,  reads  the  interface  file  and  calculates  averages  of 
each  of  the  five  variables  above  for  each  time  stored  in  the 


file  ''28:6).  SOFEPL  can  also  calculate  the  mean  and  vari¬ 
ance  of  additional  variables  which  represent  the  difference 
between  two  variables  in  he  file.  For  example,  the  differ¬ 
ence  between  a  truth  state  and  its  corresponding  filter 
state  is  an  error  state.  SOFEPL  can  calculate  the  mean  and 
variance  of  all  error  states  for  each  measurement  time. 

These  statistics  can  be  calculated  over  the  entire  Monte 
Carlo  study  or  over  a  user-defined  subset  of  runs. 

Finally,  SOFEPL  displays  the  statistics  on  computer¬ 
generated  plots.  A  number  of  different  types  of  plots  are 
available  from  SOFEPL.  Two  of  these  plots  are  of  primary 
interest  in  this  study  and  are  discussed  below. 

The  first  plot  contains  five  curves  which  display 

several  statistics  associated  with  a  single  error  state 

averaged  over  an  entire  Monte  Carlo  study.  Three  of  the 

five  curves  depict  the  actual  mean-error  of  the  filter 

state,  the  mean-error  plus  one  standard  deviation  of  the 

actual  error,  0  ,  and  the  mean-error  minus  0  .  The 

’error  error 

other  two  curves  show  the  one  standard  deviation  bounds  of 
the  error  state  as  computed  by  the  filter.  These  latter 
bounds  have  the  same  values  as  the  positive  and  negative 
square  roots  of  the  appropriate  diagonal  element  of  the 
filter  covariance  matrix  and  are  thus  labelled  as  ±  square 
roots  of  P^j(t^).  Fig.  5.1  illustrates  this  type  of  plot. 
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Several  advantages  motivate  the  use  of  this  plot.  It  en¬ 
ables  the  analyst  to  quickly  assess  overall  filter  perfor¬ 
mance  by  making  filter  divergence  or  bias  in  any  state 
immediately  apparent.  Peak  error  magnitudes,  settling  time, 
and  time-averaged  error  magnitudes  are  easily  estimated  from 
this  plot.  Additionally,  causal  relationships  are  often 
visible.  For  example,  sudden  changes  in  error  magnitude  or 
bias  might  be  associated  with  specific  events  in  the  trajec¬ 
tory,  such  as  a  fighter  roll  maneuver.  Finally,  the  inclu¬ 
sion  of  the  f i 1 te r -computed  covariance  curves  allows  the 
analyst  to  compare  the  filter's  actual  error  variance  for 
each  state  with  its  internally  calculated  representation  of 
the  error  variance  for  that  state.  This  information  pro¬ 
vides  valuable  insights  regarding  the  filter's  relative 
weighting  of  measurement  information  versus  information  from 
the  dynamics  model  (1:337-339). 

Another  valuable  SOFEPL  plot  displays  the  mean  and  one 
standard  deviation  bounds  of  a  measurement  residual.  Resid¬ 
ual  plots  provide  valuable  insights  for  the  implementation 
of  adaptive  filtering  techniques.  For  instance,  a  sudden 
increase  in  residual  values  during  a  fighter  maneuver  might 
motivate  a  shift  in  the  relative  weighting  of  measurement 
and  dynamics  information  for  the  duration  of  the  maneuver. 
Such  a  shift  could  be  achieved  by  decreasing  the  strength 
of  the  process  driving  noise  or  by  increasing  the  strength 
of  measurement  noise.  Fig.  5.2  illustrates  a  typical  SOFEPL 
residual  plot . 


Full  descriptions  of  the  two  plots  presented  above  and 


other  plots  generated  by  SOFEPL  are  available  in  Reference 
28.  Additionally,  the  interested  reader  is  referred  to 
References  13,  30,  and  31  for  a  more  detailed  discussion  of 
the  mechanics  of  SOFE  and  SOFEPL.  SOFE  and  SOFEPL  provide 
only  the  overall  Monte  Carlo  structure  for  this  study.  Com 
pleting  the  simulation  requires  an  additional  computer  pro¬ 
gram  to  serve  as  truth  model  for  the  F-4E/G  problem. 

Truth  Model 

As  explained  in  the  previous  section,  the  validity  of 
Monte  Carlo  analysis  for  any  Kalman  filter  application  is 
highly  dependent  on  the  quality'of  the  truth  model  used  as 
the  standard  of  comparison  for  filter  performance.  Ob- 


(2)  The  program  was  not  designed  for  Kalman  filter 
analysis  and  did  not  directly  calculate  relative 
target  position,  target  acceleration,  and  other 
quantities  needed  from  the  truth  model. 

(3)  The  program  only  accommodated  two  dimensional 
trajectories  for  the  attacker  and  target  aircraft, 
a  deficiency  which  led  to  observability  problems 
in  the  three  dimensional  filters. 

Modifying  the  truth  model  software  to  overcome  these  prob¬ 
lems  required  considerably  more  effort  than  originally  esti¬ 
mated.  However,  no  better  model  was  available.  This  fact 
made  program  modifications  unavoidable. 

The  overall  goal  of  the  modified  program  is  to  produce 
realistic  position,  velocity,  and  acceleration  histories  for 
typical  LRI  combat  trajectories  and  to  preserve  the  model¬ 
ling  of  dynamic  lags  in  the  APQ-120  gimbal  angle  and  angle 
rate  measurements  which  result  from  such  trajectories. 
Additionally,  the  truth  model  must  generate  ownship  attitude 
and  airspeed  information  for  the  maneuvering  F-4E/G.  Al¬ 
though  produced  by  a  two  dimensional  model,  all  variables 
must  be  oriented  in  the  three  dimensional  geometry 
of  the  Kalman  filters  under  evaluation.  Finally,  the  entire 
set  of  time-tagged  truth  states,  radar  measurements,  air¬ 
speed  values,  and  deterministic  attitude  angles  must  be 
written  into  an  unformatted  external  trajectory  file  com¬ 
patible  with  the  SOFE  simulation. 


Following  the  modification  and  verification  of  the 
truth  model  program,  three  standard  trajectories  were  de¬ 
veloped  for  Kalman  filter  evaluations.  All  three  trajec¬ 
tories  incorporate  an  F-4E/G  roll  maneuver  as  a  critical 
performance  test.  Each  trajectory  presents  a  target  located 
at  an  initial  range  of  100,000  feet  and  constrained  to  the 
same  altitude  as  the  F-4E/G.  To  remain  consistent  with  the 
LRI  target  assumptions  discussed  in  Chapter  II,  the  target 
maintains  a  non-maneuvering  flight  path  and  constant  veloc¬ 
ity  of  800  feet  per  second.  The  trajectories  vary  according 
to  the  aspect  and  velocity  which  the  target  presents  rela¬ 
tive  to  the  F-4E/G. 

The  first  scenario,  labelled  the  beam  trajectory,  simu¬ 
lates  a  westbound  target  crossing  right  to  left  in  front  of 
a  northbound  F-4E/G.  The  target's  initial  position  is 
northeast  of  the  F-4E/G  on  a  line  45  degrees  to  the  right  of 
the  attacker's  nose.  The  F-4E/G  maintains  a  constant  veloc¬ 
ity  of  800  feet  per  second.  This  is  slightly  more  than  the 
nominal  combat  airspeed  of  760  feet  per  second  listed  in 
Table  II-2,  but  is  used  in  this  simulation  to  provide  data 
for  comparison  with  other  studies  conducted  using  that  air¬ 
speed  (8,32).  After  three  seconds  of  straight  and  level 
flight,  the  F-4E/G  begins  a  one  radian  per  second  roll  to 
the  right.  At  four  seconds,  the  aircraft  stabilizes  its 
roll  angle  and  holds  a  constant  bank  turn  until  five  sec¬ 
onds.  At  that  time,  the  fighter  initiates  a  one  radian  per 
second  rollout  maneuver  which  achieves  wings  level  flight  at 


approximately  six  seconds  total  time  elapsed.  The  F-4E/G 
maintains  the  resulting  heading  and  does  not  maneuver  for 
the  remainder  of  the  20  second  trajectory.  Fig.  5.3  illus¬ 
trates  the  beam  trajectory. 

The  second,  or  tailchase,  trajectory  simulates  a  target 
fleeing  to  the  northeast  with  the  F-4E/G  maneuvering  for 
pursuit,  as  shown  in  Fig.  5.4.  Initial  target  position  is 
in  front  of  the  northbound  attacker,  but  bearing  five  de¬ 
grees  to  the  right  of  the  attacker's  nose.  The  airspeed  of 
the  F-4E/G  is  again  constant,  but  set  at  850  feet  per  second 
in  this  trajectory  to  permit  the  attacker  to  achieve  a 
slight  closure  rate  on  the  target.  This  prevents  the  range 
from  exceeding  LRI  limits.  The  F-4E/G  remains  wings  level 
northbound  for  five  seconds  before  initiating  a  one  radian 
per  second  right  roll  and  stabilizing  in  a  constant  bank 
right  turn  one  second  later.  This  turn  is  held  for  12  more 
seconds  before  the  attacker  rolls  out  to  resume  wings  level 
flight  in  a  position  almost  directly  astern  the  target  and 
at  a  heading  which  diverges  slightly  to  the  left  (north)  of 
the  target's  heading.  This  flight  path  is  held  for  the 
remainder  of  the  20  second  trajectory  time. 

The  final  test  trajectory  simulates  a  head-on  attack 
profile.  The  target  approaches  from  the  northeast,  heading 
185  degrees.  The  northbound  F-4E/G  maintains  a  constant  800 
feet  per  second  airspeed  and  duplicates  the  maneuver  de¬ 
scribed  in  the  beam  trajectory;  a  right  roll  at  three  sec¬ 
onds,  constant  bank  at  four  seconds,  and  rollout  to  wings 
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Tail  chase  Trajectory 


Le.t.'I  jc  five  seconds.  This  maneuver  places  the  target  and 
i :  t  j  c  K  e  t  ' j  n  reciprocal  headings  with  a  5 1) 0  foot  lateral 
i  L  :  s  i' t  as  shown  in  F  ’  a  ,  5.5.  Both  aircraft  maintain  h  e  a  d  i  n  a 
!  ir  ine  Is  remainina  seconds  of  this  test. 

!'  1;  e  a  e  n  e  r  a  l  1  «,>  n  of  the  above  truth  state  trajectories 
;ri;Laili?s  the  calculation  of  corresponding  time  histories  for 
ridar  measurements,  o unship  attitude  angles,  and  owns  hip 
iirspeed.  Tiiese  trajectories  complete  the  software  tools 
needed  to  conduct  a  Monte  Carlo  analysis  of  the  CV  and  CM. A 
tillers.  Ttie  next  section  discusses  general  procedures  for 
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General  Evaluation  Procedures 


The  preceding  sections  of  this  chapter  explain  the 
theoretical  basis  of  Monte  Carlo  analysis  and  describe  the 
structure  of  an  actual  Monte  Carlo  simulation  program 
designed  to  evaluate  the  extended  Kalman  filters  developed 
in  this  study.  This  section  discusses  methods  of  actually 
using  the  above  simulation  program  to  perform  the  desired 
evaluation.  All  of  the  methods  presented  are  iterative 
procedures  which,  in  their  most  fundamental  sense,  have  the 
following  form: 

(1)  Perform  a  Monte  Carlo  study  on  a  candidate  filter 
to  obtain  a  statistical  characterization  of  filter 
performance. 

(2)  Critically  examine  the  filter  performance  from  the 
previous  step  by  comparing  it  to  (a)  the  perfor¬ 
mance  of  other  candidate  filters  operating  under 
identical  conditions,  (b)  performance  of  the  same 
filter  with  varying  noise  strengths,  time  constants 
or  other  variable  parameters,  and  (c)  predetermined 
standards  of  filter  performance. 

( ’)  Using  insights  gained  from  the  previous  step, 

adjust  the  values  of  filter  parameters  or  change 
elements  of  the  filter  structure  with  the  goal  of 
improving  filter  performance. 

(4)  Repeat  the  above  steps  until  satisfactory  perfor¬ 
mance  is  obtained. 


In  order  to  perform  meaningful  comparisons  of  the  type 
suggested  in  Step  2  of  the  above  procedure,  some  metric  or 
yardstick  of  filter  performance  is  necessary.  For  this 
study,  the  following  three  figures  of  merit  are  used: 

(1)  peak  error  magnitude  for  each  position  and  veloc¬ 
ity  state, 

(2)  settling  time,  defined  as  the  time  between  the 
onset  of  an  F-4E/G  maneuver  and  the  time  filter 
error  states  stabilize  within  five  per  cent  of 
peak  values, 

and  , 

(3)  time-averaged  error  magnitude  for  each  position 
and  velocity  state  estimated  over  the  length  of  a 
given  trajectory.  This  value  has  little  meaning 
except  in  relation  to  a  specific  trajectory  and  is 
used  only  to  compare  performance  among  filter  ver¬ 
sions  operating  over  identical  trajectories. 

The  above  figures  of  merit  are  chosen  because  they 
characterize  filter  performance  in  terms  consistent  with  the 
performance  criteria  for  an  adequate  LRI  filter  as  identi¬ 
fied  in  Chapter  I.  Additionally,  these  figures  of  merit  are 
also  easily  obtained  from  available  SOFEPL  plots.  Accelera¬ 
tion  states  are  not  included  in  the  figures  of  merit  because 
they  provide  very  little  information  in  the  benign  trajec¬ 
tories  associated  with  LRI  targets,  and  are  not  available  in 
the  case  of  the  CV  model.  Acceleration  states  are  examined. 


however,  for  evidence  of  filter  divergence  and/or  observa¬ 
bility  problems,  as  discussed  in  Chapter  VI. 

The  above  figures  of  merit  can  be  used  in  several  ways 
to  provide  insights  for  improving  filter  performance  over  a 
succession  of  Monte  Carlo  studies.  During  the  initial  phase 
of  filter  evaluation,  for  instance,  the  filter  designer 
seeks  to  validate  filter  models  and  the  software  which 
implements  these  models.  One  way  of  doing  this  is  to  use 
comparisons  of  figures  of  merit  to  establish  consistency  in 
the  baseline  performance  of  the  filters  under  evaluation. 

In  this  sense,  consistency  implies  that  several  similarly 
structured  filters  or  filter  variations  are  likely  to  ex¬ 
hibit  the  same  general  performance  patterns  when  tested  on 
the  same  trajectory.  When  comparing  filter  performance 
patterns,  the  analyst  expects  variations  in  peak  error  and 
other  figures  of  merit.  If  the  figures  of  merit  vary  in  a 
consistent  manner,  overall  confidence  in  the  filter  simula¬ 
tion  is  enhanced.  However,  inconsistencies  in  filter  per¬ 
formance  often  point  to  serious  modelling  errors  or  mistakes 
in  computer  code. 

As  an  example,  consider  two  Monte  Carlo  studies  of  the 
GMA  filter  operating  in  the  beam  trajectory  scenario.  The 
first  study  incorporates  only  four  of  the  available  measure¬ 
ments;  range,  range  rate,  and  the  sines  of  the  azimuth  and 
elevation  angles.  The  second  study  differs  only  in  that  it 
incorporates  all  six  measurements  into  the  state  estimates. 
Comparing  the  results  of  these  two  studies,  the  analyst 


finds  a  close  resemblence  in  the  general  patterns  of  filter 
performance  and  perhaps  notes  a  modest  decrease  in  peak 
error  magnitude  for  the  second  filter,  indicating  that  the 
additional  measurements  improve  filter  performance  in  that 
respect.  A  series  of  consistent  results  of  this  nature 
lends  confidence  to  the  GMA  model  and  its  computer  implemen¬ 
tation.  However,  if  results  of  the  two  studies  show  stable 
convergent  performance  in  the  first  case,  and  several  rap¬ 
idly  diverging  states  in  the  second  case,  then  a  problem 
exists  which  the  analyst  must  investigate.  This  problem  is 
likely  related  to  the  angle  rate  measurements,  but  might  be 
in  the  form  of  a  mismodelled  measurement,  a  sign  error  in  a 
filter  measurement  equation,  an  error  in  the  truth  model,  a 
typographical  error  in  the  computer  code,  or  some  combina¬ 
tion  of  these.  Examples  of  each  of  these  errors  were  actu¬ 
ally  uncovered  during  preliminary  Monte  Carlo  studies  for 
the  F-4E/G  problem. 

Other  examples  of  filter  evaluation  procedures  include 
filter  "tuning"  and  the  determination  of  an  error  budget. 
Brief  descriptions  of  these  activities  are  given  below. 

More  detailed  examples  of  their  use  and  the  results  of  their 
application  in  this  study  are  discussed  in  Chapter  VI. 

Filter  "tuning"  describes  the  process  of  iteratively 
adjusting  the  value  of  a  single  filter  parameter  based  on 
results  of  successive  Monte  Carlo  studies.  Commonly,  the 
parameter  chosen  for  adjustment  is  Q(t),  the  strength  of 
dynamics  driving  noise,  since  its  value  is  seldom  known  with 


great  certainty.  By  repeating  Monte  Carlo  studies  and  using 
different  values  of  Q(t),  the  filter  analyst  can  find  driv¬ 
ing  noise  strenths  which  result  in  lower  peak  errors, 
shorter  settling  times,  or  other  desirable  qualities  as 
measured  by  figures  of  merit.  Similarly,  a  series  of  Monte 
Carlo  studies  can  be  conducted  by  leaving  Q(t)  constant  and 
varying  another  parameter,  such  as  measurement  noise 
strengths  or,  in  the  case  of  the  GMA  model,  the  target 
acceleration  time  constant,  T. 

An  error  budget  is  a  technique  used  to  determine  the 
effect  of  separate  sources  of  error  on  filter  performance 
(1:339).  If  each  known  source  of  error  is  modelled  sepa¬ 
rately  in  the  truth  model,  it  is  posssible  to  turn  on  the 
error  sources  one  at  a  time  for  a  series  of  Monte  Carlo 
studies.  Using  the  figures  of  merit  as  metrics,  the  analyst 
compares  the  results  of  these  Monter  Carlo  studies  and 
assesses  the  relative  contribution  of  each  separate  error 
source  to  the  overall  filter  error.  By  identifying  the 
largest  contributors  to  filter  error,  the  analyst  determines 
areas  in  which  adjustments  or  changes  in  filter  design  are 
most  likely  to  yield  significant  performance  improvements. 

Conceptually,  the  above  filter  evaluation  procedures 
can  be  separated  into  distinct  activities  and  discussed  in  a 
logical  sequence.  In  practice,  however,  these  activities 
are  strongly  interdependent  and  are  frequently  accomplished 
iteratively  or  in  parallel.  For  example,  an  error  budget  is 
typically  constructed  using  a  filter  which  has  already  been 


validated  and  tuned.  However,  the  results  of  the  error 


budget  may  mandate  filter  changes  which  require  re-valida¬ 
tion  and  retuning  of  the  modified  filter.  In  turn,  another 
error  budget  may  be  justified.  The  intricacies  of  the 
filter  evaluation  process  vary  from  application  to  applica¬ 
tion  and  are  best  understood  in  the  context  of  a  specific 
problem.  The  next  chapter  discusses  evaluation  procedures 
as  applied  to  the  proposed  F-4E/G  LRI  filters  and  presents 
the  results  from  each  stage  of  the  filter  evaluation 
process . 
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Introduction 

This  chapter  presents  the  results  of  Monte  Carlo  stud¬ 
ies  conducted  to  evaluate  the  performance  of  the  CV  and  GMA 
filters  for  the  three  test  trajectories  defined  in  the  prev 
ious  chapter.  Initial  Monte  Carlo  studies,  which  are  used 
to  identify  and  correct  flaws  in  the  simulation  software, 
are  not  included  here  because  they  do  not  provide  an  accu¬ 
rate  portrayal  of  actual  filter  performance.  Later  studies 
conducted  after  the  establishment  of  a  credible  simulation 
package,  provide  a  better  representation  of  filter  perfor¬ 
mance  and  are  used  for  filter  tuning  and  other  evaluation 
procedures.  These  studies  are  presented  sequentially  as 
f  o 1 lows : 

(1)  Baseline  filter  comparisons  and  validation 

(2)  Filter  tuning  using  driving  noise 

(3)  Error  budget  studies 

(4)  Filter  tuning  using  measurement  noise 

(5)  Filter  adaptation  to  F-4E/G  roll  rate 

(6)  Filter  comparisons  on  three  trajectories 

As  emphasized  in  the  previous  chapter,  filter  evalua¬ 
tion  procedures  are  highly  interrelated.  Most  of  the  above 
studies  are  repeated  several  times  in  response  to  the  out¬ 
comes  of  related  studies.  In  most  cases,  the  results  tabu¬ 
lated  in  this  chapter  are  the  culmination  of  repetitive 


testing.  Furthermore,  the  order  of  the  above  activities 
does  not  imply  their  exact  order  of  accomplishment. 


Baseline  Filter  Comparisons  and  Validation 

One  of  the  early  goals  of  the  filter  evaluation  process 
is  to  obtain  side-by-side  comparisons  of  the  CV  and  GMA 
filters  operating  under  comparable  conditions.  To  attain 
this  goal,  a  five-run  Monte  Carlo  study  is  conducted  for 
each  filter  using  the  beam  trajectory.  The  beam  trajectory 
considered  the  worst  case  trajectory  because  the  los 
rates  in  this  trajectory  caused  lags  in  the  antenna  angle 
and  angle  rate  measurements  which  produced  larger  errors 
than  in  the  other  two  trajectories.  The  beam  trajectory  is 
also  the  primary  test  case  used  for  filter  tuning  and  other 
evaluation  procedures.  The  baseline  >  >mparison  studies  use 
the  initial  strengths  of  driving  noise  derived  in  Chapter  II 
for  the  CV  and  GMA  filters.  Initial  measurement  noise 
strengths  were  provided  by  F-4E/G  software  engineers  at  Hill 
AFB,  Utah,  as  discussed  in  Chapter  III. 

Although  0.04  second  is  the  original  update  interval 
specified  in  Chapter  I,  allowance  is  made  for  longer  update 
intervals.  The  maximum  allowable  interval  is  0.1  second,  a 
constraint  imposed  by  the  10  Hz  rate  required  for  cockpit 
instrumentation  updates.  In  terms  of  filter  accuracy,  the 
0.1  second  update  rate  is  considered  the  most  severe  condi¬ 
tion  since  it  incorporates  the  least  measurement  information 
and  allows  propagation  error  to  grow  uncorrected  for  the 


longest  period  of  time.  For  this  reason,  0.1  second  is  the 
measurement  interval  used  for  all  filter  evaluation  studies. 
Theoretically,  an  actual  filter  implementation  which  uses  a 
faster  update  rate  could  achieve  improved  performance  after 
recalculation  of  ,  as  given  by  (4-16). 

The  results  of  the  baseline  filter  comparisons  are 
illustrated  in  Appendix  E.  In  studies  using  all  six  meas¬ 
urements,  both  filters  performed  well  until  the  initiation 
of  the  F-4E/G  roll  maneuver,  when  both  filters  produced 
large  errors  in  position  and  velocity  states.  The  GMA 
filter  converged  to  pre-maneuver  accuracy  in  all  states 
within  six  seconds  after  the  end  of  the  maneuver.  The  CV 
filter  produced  errors  an  order  of  magnitude  greater  than 
the  GMA  filter  and  failed  to  converge  in  all  states  after 
the  maneuver  ended.  These  results  are  illustrated  by  SOFEPL 
plots  in  Appendix  E,  Figs.  E.l  through  E.15.  Table  VI-1 
summarizes  the  results. 

A  second  set  of  studies  compared  CV  and  GMA  filter 
performance  using  only  four  measurements;  range,  range  rate, 
and  the  two  gimbal  angles.  Again,  both  filters  committed 
unacceptably  large  position  and  velocity  errors  during  own- 
ship  maneuvers,  but  in  this  case,  both  filters  exhibited 
errors  of  similar  magnitudes  and  displayed  settling  times  of 
7-12  seconds  in  all  states.  Table  VI-2  summarizes  these 
studies.  Full  performance  plots  are  shown  in  Figs.  E.16 
through  E.20.  Since  the  peak  total  velocity  errors  commit¬ 
ted  by  both  filter  in  these  initial  studies  are  two  orders 


State 

Peak 

CV 

Error 

Magnitude 

GMA 

Settling 

CV 

time 

GMA 

Position 

North 

65000 

ft 

980 

ft 

none 

10.5 

sec 

East 

61000 

ft 

980 

ft 

none 

12.0 

sec 

Down 

35000 

ft 

3500 

ft 

none 

11.6 

sec 

Velocity 

North 

41000 

f  ps 

590 

f  ps 

10.0  sec 

10.5 

sec 

East 

35000 

f  ps 

710 

f  ps 

9 . 0  sec 

10.3 

sec 

Down 

27000 

f  ps 

2800 

f  ps 

9 . 0  sec 

11.0 

sec 

TABLE  VI-2 


Baseline  Comparison  of  CV  and  GMA  Filters  on  Beam  Trajectory 

(4  Measurements) 


State 

Peak 

Error 

Magnitude 

Settling 

time 

CV 

GMA 

CV 

GMA 

Position 

North 

1250 

ft 

1125  ft 

7.3  sec 

10.2  sec 

East 

1270 

ft 

1100  ft 

11.5  sec 

10.0  sec 

Down 

4400 

ft 

3300  ft 

7 . 3  sec 

11.5  sec 

Velocity 

North 

1090 

f  ps 

475 

f  ps 

8.8 

sec 

10.2 

sec 

East 

1250 

f  ps 

780 

f  ps 

9.0 

sec 

11.9 

sec 

Down 

3900 

f  ps 

2520 

f  ps 

9.6 

sec 

11.0 

sec 

of  magnitude  higher  than  desired,  filter  tuning  is  required. 
An  examination  of  the  above  data  reveals  that  rhe  GMA  filter 
with  four  measurements  committed  the  lowest  peak  total 


velocity  error,  2680  fps.  For  this  reason  most  filter 
evaluation  procedures  after  this  point  use  only  four  meas¬ 
urements;  range,  range  rate,  and  gimbal  angle  sines. 

Comparison  of  CV  and  GMA  filter  performance  using  four 
measurements  supports  the  expected  result  that  the  GMA  fil¬ 
ter  provides  improved  estimation  over  the  CV  model.  Peak 
position  errors  were  10-25  %  lower  and  peak  velocity  errors 
were  35-56  %  lower  for  the  GMA  filter  as  compared  to  the  CV 
filter.  The  relatively  large  errors  in  the  vertical  channel 
are  attributed  to  a  lack  of  observability  in  that  dimension, 
a  consequence  of  the  two-dimensional  truth  model.  Settling 
times  were  slightly  longer  for  the  GMA  model,  but  remained 
less  than  12  seconds.  Time-averaged  error  is  not  used  for 
this  portion  of  the  evaluation  because  it  provides  little 
useful  information  when  peak  error  magnitudes  are  very  large. 

Filter  Tuning  Using  Driving  Noise 

This  series  of  Monte  Carlo  studies  sought  to  improve 
baseline  performance  of  the  GMA  filter  by  adjusting  the 
strength  of  dynamic  driving  noise.  Several  studies  were 
conducted  using  the  GMA  filter  with  the  four  measurements 
discussed  in  the  previous  section.  These  studies  indicate 
that  filter  performance  improves  as  Q(t)  is  decreased  from 
the  baseline  value.  Table  VI-3  summarizes  a  series  of 
simulations  in  which  Q(t)  is  lowered  by  factors  of  10  and 
compared  with  baseline  filter  performance.  Appendix  F  con¬ 
tains  performance  plots  for  these  studies. 

106 


TABLE  VI-3 


Effect  of  Changing  Driving  Noise  Strength  on  GMA  Filter 

(4  Measurements) 


Q(t)  : 

(ft**  / sec^  ) 

5000 

1000 

100 

10 

Position : 

Peak  Error 

Magnitude  ( f t ) /Se t t 1 i ng  time  (sec) 

North 

East 

Down 

1  125/10.2 

1  100/10.0 
3300/  1  1  . 5 

1000/12.0 

990/12.5 

2800/12.5 

970/13.5 

980/15.0 

2450/13.0 

9  50/* 
900/17.0 
2400/* 

Velocity  ; 

Peak  Error 

Magnitude  (fp 

s )/Settling 

time  ( sec  ) 

North 

East 

Down 

475/10.2 
780/  1  1.9 
2520/ 11.0 

360/12.3 

365/12.5 

1450/12.0 

265/17.0 

250/16.5 

770/13.5 

250/* 

240/* 

710/* 

Greater  than  20  sec 


Two  trends  are  apparent  in  Table  VI-3.  As  Q(t)  is 
decreased,  peak  error  magnitudes  decrease  and  settling  times 
increase.  A  value  of  100  ft‘‘/sec®  is  chosen  as  a  good 
value  of  Q( t ) ,  since  it  results  in  reasonable  values  for 
both  peak  error  and  settling  time.  This  value  is  used  in 
the  remaining  filter  evaluations. 

Error  Budget  Studies 

Although  the  truth  model  used  for  filter  evaluation 
does  not  model  individual  error  sources  in  the  dynamic  sys¬ 
tem,  it  is  possible  to  investigate  the  effect  of  different 
measurement  noise  sources  by  running  the  filter  with  differ¬ 
ent  combinations  of  one,  two,  or  three  measurements.  One 


complicating  factor  that  results  from  testing  the  filter 
with  fewer  than  four  measurements  is  that  the  system  is  no 


longer  completely  observable.  In  Monte  Carlo  studies  using 
the  GMA  filter  with  three  or  fewer  measurements,  it  is 
difficult  to  separate  the  effect  of  unobservable  states 
from  other  effects,  such  as  the  effect  of  an  improperly 
modelled  measurement.  However,  useful  insights  can  be  ob¬ 
tained  by  this  technique  and  it  is  used  extensively  in  this 
study  . 

Three  general  performance  patterns  result  from  Monte 
Carlo  studies  using  less  than  a  full  complement  of  four 
measurements.  Those  studies  using  only  range,  range  rate, 
or  both  are  completely  insensitive  to  ownship  maneuver. 

Mean  errors  remain  close  to  zero  in  all  states,  but  filter- 
computed  covariances  increase  without  bound.  In  Appendix  G, 
Figs.  G.l  through  G.9  illustrate  this  performance  pattern. 

A  different  pattern  results  from  using  only  gimbal  angle 
measurements.  The  filter  error  is  extremely  sensitive  to 
ownship  maneuvers.  As  seen  in  Figs.  G.IO  through  G.18,  the 
overall  performance  of  filters  using  only  angle  measurements 
resembles  the  general  pattern  of  performance  with  four  meas¬ 
urements,  but  with  larger  peak  errors  and  longer  settling 
times  in  position  and  velocity  states.  A  third  general 
pattern  of  performance  emerges  from  studies  using  only  angle 
rate  measurements.  In  these,  several  GMA  states  diverge 
rapidly  after  ownship  maneuver,  as  illustrated  in  Figs.  G.19 
through  G.27.  This  result  gives  additional  impetus  to  con¬ 
tinue  testing  without  using  angle  rate  measurements. 


In  order  to  further  investigate  the  effect  of  dynamic 
lags  in  the  girabal  angle  measurements,  the  truth  model  is 


modified  to  remove  these  lags.  This  allows  a  critical 
source  of  error  to  be  turned  on  and  off  in  the  truth  model 
A  Monte  Carlo  study  is  performed  using  the  GMA  filter  with 
range,  range  rate,  and  gimbal  angle  measurements  with  the 
lags  removed.  The  results  show  peak  error  improvements  of 
an  order  of  magnitude  or  greater  in  all  states,  as  can  be 
seen  by  comparing  Figs.  G.28  through  G.36  with  Figs.  F.IO 
through  F.18.  Table  VI-4  summarizes  this  comparison. 

TABLE  VI-4 

Effect  of  Removing  Dynamic  Lags  from  Angle  Measurements 
(GMA  model,  4  measurements,  Q  =  100) 


With  Lags  Without  Lags 


Position:  Peak  Error  Magnitude  ( f t ) /Se t t 1 i ng  time  (sec) 

North  970/13.5  30/0.0 
East  980/15.0  30/0.0 
Down  2450/13.0  35/0.0 


Velocity:  Peak  Error  Magnitude  (fps)/ Settling  time  (sec) 

North  265/17.0  8/0.0 
East  250/16.5  8/0.0 
Down  770/13.5  14/0.0 


The  results  displayed  in  Table  VI-4  suggest  that  the 
dynamic  lags  in  the  angle  measurements  contribute  signifi¬ 
cantly  to  GMA  filter  error.  One  way  to  improve  filter 
performance  is  to  model  the  lags  in  the  filter  dynamics 


equations.  Another  way  is  to  adjust  the  value  of  measure¬ 
ment  noise  associated  with  angle  measurements,  as  discussed 
in  the  next  section. 

Filter  Tuning  Using  Measurement  Noise 

The  results  of  the  error  budget  studies  discussed  in 
the  previous  section  suggest  that  the  angle  and  angle  rate 
measurements  make  the  largest  contributions  to  overall  fil¬ 
ter  error.  The  lags  in  the  antenna  angle  measurements  are 
isolated  and  shown  specifically  to  be  an  important  source  of 
filter  error.  These  results  motivate  the  next  series  of 
''lonte  Carlo  studies,  in  which  the  strengths  of  the  measure¬ 
ment  noise  associated  with  the  angle  and  angle  rate  measure¬ 
ments  are  increased.  The  first  study,  illustrated  in  Appen¬ 
dix  li ,  Figs.  H.l  through  H.13,  evaluates  the  GMA  filter 
using  four  measurements  and  a  driving  noise  strength  of  100 
ft  “/sec  This  allows  direct  comparison  with  the  filter 
performance  study  shown  in  Figs.  F.IO  through  F.18,  in  which 
the  original  angle  measurement  noise  strength  of  1.1  x  10  ^ 
is  used  for  both  angle  measurements.  In  the  new  study, 

measurement  noise  strengths  for  the  two  angle  measurements 

-4 

are  increased  to  J.O  x  10  ,  as  discussed  below.  The  in¬ 

crease  in  measurement  noise  yields  improvements  of  76.4  to 
MO. 8  %  in  the  position  states  and  improvements  of  87.5  to 
')  1  . 0  T  in  the  velocity  states.  Settling  times  increase  2  to 
4  seconds  in  position  states,  but  decrease  a  similar  amount 


TABLE  VI-5 


Effect  of  Increasing  Measurement  Noise  Strengths  Associated 
With  Antenna  Angle  Measurements 

(GMA  model,  4  measurements,  Q  =  100) 


Angle  Measurement 
Noise  Strength: 

1.1  X  10  ^ 

2.0  X  10"^ 

Position: 

Peak 

Error  Magnitude  ( f t ) /Se t t 1 i ng  time  (sec) 

North 

970/13.5 

230/15.0 

East 

980/15.0 

230/19.0 

Down 

2450/13.0 

470/13.0 

Velocity  : 

Peak 

Error  Magnitude  (fps 

0/Settling  time  (sec) 

North 

265/17.0 

33/14.3 

East 

250/16.5 

29/13.0 

Down 

770/13.5 

62/13.0 

The  value  of  angle  measurement  noise  strength  in  the 
right  column  Table  VI-5  was  selected  after  several  earlier 
Monte  Carlo  studies  indicated  that  higher  values  did  not 
produce  significant  improvements  in  filter  performance.  The 
lowest  possible  value  is  desirable  because  the  filter  essen¬ 
tially  ignores  a  measurement  if  its  associated  measurement 
noise  strength  is  too  large.  The  value  of  2.0  x  10~^  also 
agreed  closely  with  a  computer  experiment  in  which  the  truth 
model  was  adapted  to  calculate  the  variance  of  the  sine  of 
the  angle  errors  during  an  F-4E/G  roll  maneuver. 

Another  group  of  Monte  Carlo  studies  investigates  the 
performance  of  the  GMA  filter  with  six  measurements  after 
increasing  the  noise  strength  of  the  angle  rate  measure- 


merits.  The  results  of  this  study,  shown  in  Figs,  H.14 
through  H.28,  indicate  that  no  performance  advantage  is 
gained  by  incorporating  the  angle  rate  measurements.  Peak 
error  statistics  duplicate  the  right  column  of  Table  VI-5 
except  for  an  8-10  %  reduction  in  the  vertical  channel 
position  and  velocity  states.  Settling  times  increase  in 
the  six  measurement  simulation,  and  exceed  the  20  second 
simulation  time  for  the  vertical  channel  states.  For  these 
reasons,  no  further  evaluations  are  conducted  using  six 
measurements.  Other  techniques  for  improving  GMA  filter 
performance  are  discussed  in  the  next  section. 

Filter  Adaptation  to  F -4E/G  Roll  Rate 

Adaptive  filtering  techniques  are  investigated  in  an 
attempt  to  further  improve  filter  performance.  Although 
adaptive  filtering  is  typically  based  on  residual  monitoring 
schemes  (2:69),  a  more  direct  method  is  chosen  for  this 
study.  Since  the  most  serious  filter  errors  occur  in  re¬ 
sponse  to  F-4E/G  maneuvers,  a  simple  attitude  rate  detector 
can  be  used  to  trigger  filter  adaptations.  Two  studies  are 
conducted  on  the  GMA  model  using  this  approach.  In  the 
first,  a  roll  detector  increases  the  value  of  gimbal  angle 
measurement  noise  whenever  ownship  roll  rates  exceed  a 
threshold  of  0.05  radians  per  second.  In  the  second  study, 
the  value  of  Q(t)  is  decreased  to  1  ft  "*/sec  ^  when  roll  rate 
exceeds  the  threshold.  Both  studies  produced  disappointing 
results.  The  adaptive  Q(t)  filter  shows  little  change  in 


performance,  while  the  adaptive  measurement  noise  filter 
commits  peak  errors  four  times  larger  than  the  same  filter 
without  adaptation.  In  both  cases,  perfotmance  plots  indi¬ 
cate  that  adaptations  occur  at  the  desired  times  but  do  not 
produce  the  desired  effects.  Appendix  I  contains  perform¬ 
ance  plots  for  these  studies. 

Filter  Comparisons  on  Three  Tra  jectories 


As  a  final  demonstration  of  the  feasibility  of  a  Kalman 
filter  for  the  F-4E/G  LRI  function,  the  tuned  GMA  filter  is 
compared  in  all  three  trajectories  discussed  in  Chapter  V. 
The  comparisons  involve  20-run  Monte  Carlo  simulations  of 
the  GMA  filter  using  four  measurements.  The  results  indi¬ 
cate  similar  performance  for  the  beam  and  tailchase  trajec¬ 
tories.  Performance  improves  significantly  in  the  head-on 
attack  profile.  Appendix  J  contains  performance  plots  of 
these  studies,  which  are  summarized  in  Table  VI-6. 

With  the  exception  of  the  head-on  attack  profile,  the 
final  evaluations  continue  to  display  vertical  channel 
errors  approximately  twice  as  large  as  errors  in  the  other 
channels.  As  stated  earlier,  this  is  attributed  to  a  lack 
of  observability  resulting  from  the  two  dimensional  truth 
model  trajectories.  An  attempt  was  made  to  evaluate  GMA 
filter  performance  while  masking  the  lack  of  observability 
in  the  vertical  channel.  This  was  done  by  repeating  the 
20-run  Monte  Carlo  studies  of  each  trajectory  discussed 
above  and  artificially  suppressing  the  errors  in  all  three 


TABLE  VI-6 


GMA  Filter  Performance  on  Three  Trajectories 
(4  measurements) 


BEAM 

TAILCHASE 

HEAD-ON 

Position  : 

Peak/Time-averaged  Error 

Magnitudes  ( f  t ) 

North 

175/75 

190/50 

10/9 

East 

175/78 

203/60 

65/30 

Down 

345/90 

500/100 

50/25 

Velocity : 

Peak /Time- 

averaged  Error 

Magnitudes  (fps) 

North 

24/8 

23/5 

2/0 

East 

23/10 

25/5 

8/2 

Down 

45/4 

58/12 

7/3 

vertical  channel  states.  The  results  of  this  ad  hoc  proce¬ 
dure  were  inconclusive.  North  and  east  channels  produced 
approximately  the  same  peak  errors  and  settling  times  ob¬ 
served  without  vertical  channel  supression.  In  some  states, 
peak  errors  were  2-5  Z  larger  than  in  full  state  perfor¬ 
mance.  Appendix  K  contains  plots  of  these  studies. 

An  analysis  of  the  target  velocity  errors  in  Table  VI-6 
shows  peak  total  velocity  error  values  of  55.9  fps  for  the 
beam  trajectory,  67.2  fps  for  the  tailchase,  and  10.8  fps 
for  the  head-on  attack.  If  these  values  are  considered  to 
be  three-sigma  values  of  total  velocity  error,  the  one 
standard  deviation  values  of  velocity  error  are  18.6,  22.4, 
and  3.6  fps,  respectively.  In  each  case  this  is  lower  than 
the  30  fps  benchmark  established  in  Chapter  I,  and  suggests 


TABLE  VI-7 


Comparison  of  Tuned  CV  and  GMA  Filters 
(4  measurements,  20  runs.  Beam  Trajectory) 


CV 

GMA 

Position : 

Peak /Time-averaged 

Error  Magnitudes 

(ft) 

North 

180/80 

175/75 

East 

180/90 

175/78 

Down 

350/50 

345/90 

Velocity: 

Peak /Time-averaged 

Error  Magnitudes 

(  f  ps ) 

North 

28/10 

24/8 

East 

26/11 

23/10 

Down 

50/20 

45/4 

that  the  GMA  filter  with  four  measurements  can  exceed 
Wiener-Hopf  filter  performance. 

A  final  validation  of  the  CV  filter  is  also  accom¬ 
plished  using  four  measurements  and  noise  values  obtained 
from  the  previous  filter  evaluations.  A  20-run  Monte  Carlo 
analysis  of  CV  filter  performance  on  the  beam  trajectory  is 
illustrated  in  Appendix  D,  Figs.  D.l  through  D.6  and  summar¬ 
ized  in  Table  VI-7.  The  overall  standard  deviation  of 
velocity  error  for  the  CV  filter  is  21.0  fps  in  this  study, 
which  establishes  this  filter  as  a  viable  alternative  for 


use  in  the  F-4E/G  LRI  function. 


VII,  Cone lusions  and  Recommendations 

As  stated  in  the  introductory  chapter,  the  purpose  of 
this  study  is  to  determine  the  feasibility  of  improving  the 
LRI  function  of  the  F-4E/G  fire  control  system  by  replacing 
the  existing  Wiener-Hop  filter  with  a  Kalman  filter.  To 
warrant  implementation,  the  proposed  Kalman  filter  must 
produce  target  velocity  estimates  which  meet  LRI  accuracy 
criteria  regardless  of  F-4E/G  combat  maneuvering.  Addition¬ 
ally,  the  Kalman  filter  must  operate  under  the  severe  con¬ 
straints  imposed  by  limited  memory,  slow  speed,  short  word- 
length,  and  noise-corrupted  measurements  available  in  the 
existing  F-4E/G  fire  control  computer.  Two  candidate  fil¬ 
ters  are  proposed,  one  based  on  a  constant  velocity  dynamics 
model,  and  the  other  based  on  a  Gauss-Markov  target  acceler¬ 
ation  model.  Both  models  are  implemented  as  extended  Kalman 
filters  in  a  computer  simulation  and  evaluated  using  Monte 
Carlo  analysis.  The  results  of  the  simulation  studies  pro¬ 
vide  the  basis  for  the  conclusions  and  recommendations  pre¬ 
sented  in  this  chapter. 

Conclusions 

A  9-state  reduced  order  Kalman  filter  algorithm  can  be 
implemented  in  the  constrained  memory  of  the  LRU-1  fire 
control  computer  and  can  operate  within  the  allowable  meas¬ 
urement  interval,  providing  that  the  requirement  for  on-line 


numerical  integration  is  avoided.  This  requirement  is  met 
by  preserving  linear,  time-invariant  dynamics  in  the  filter 
model.  Both  the  proposed  CV  and  GMA  models  are  satisfactory 
in  this  respect,  since  inherent  nonlinearities  in  these 
models  are  confined  to  the  measurement  equations.  The 
results  obtained  from  the  performance  evaluations  of  these 
filters  support  the  conclusion  that  the  CV  and  GMA  filters 
are  feasible  replacements  for  the  existing  Wiener-Hopf  fil¬ 
ter  in  the  F-4E/G  OFF.  The  GMA  filter  is  emphasized  for  two 
reasons.  First,  its  performance  exceeds  that  of  the  CV 
filter  when  evaluated  in  LRI  simulations.  Second,  varia¬ 
tions  of  the  GMA  filter  can  be  adapted  to  successfully 
perform  estimation  for  the  F-4E/G  ACM  function  as  well  as 
for  the  LRI  function  (32).  A  single  filter  capable  of 
satisfying  the  requirements  of  both  the  LRI  and  ACM  func¬ 
tions  is  strongly  motivated  by  the  speed  and  memory  limita¬ 
tions  of  the  LRU-1  computer. 

Although  the  performance  results  of  the  tuned  GMA  fil¬ 
ter  indicate  that  its  velocity  error  has  a  standard  devia¬ 
tion  well  below  that  of  the  Wiener-Hopf  benchmark  value  of 
30  fps  on  all  trajectories  tested,  there  are  a  number  of 
unresolved  issues  which  impact  the  ability  of  the  GMA  filter 
to  achieve  similar  performance  in  an  operational  setting. 
These  issues  include  the  failure  of  the  angle  rate  measure¬ 
ments  to  improve  filter  performance,  the  modelling  of  the 
dynamic  lags  in  the  sine  angle  measurements,  and  the  pres¬ 
ence  of  large  vertical  channel  errors  attributed  to  lack  of 
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observability  in  the  simulated  trajectories.  The  need  to 
address  these  issues  before  attempting  any  operational  fil¬ 
ter  implementation  motivates  many  of  the  recommendations  in 
the  next  section. 

Recommendations 

The  conclusions  reached  in  this  study  lead  to  recommen¬ 
dations  in  two  general  areas,  advanced  testing  and  opera¬ 
tional  implementation.  Recommendations  concerning  advanced 
testing  are  motivated  by  the  existence  of  the  unresolved 
issues  discussed  above  and  the  need  to  critically  examine 
their  impact  on  an  operational  filter.  Attention  should  be 
focused  on  aspects  of  both  the  filter  and  truth  models.  One 
primary  recommendation  is  to  expand  the  truth  model  to 
accommodate  three  dimensional  aircraft  trajectories.  This 
would  help  resolve  the  question  of  observability  in  the 
vertical  channel.  Another  valuable  improvement  is  to  model 
individual  sources  of  error  separately  in  the  truth  model. 
Examples  of  errors  that  need  to  be  modelled  are  the  glint 
noise  error  in  radar  measurements,  the  errors  associated 
with  INS  attitude  angles,  and  errors  in  the  ownship  velocity 
vector.  These  error  models  complement  the  previously  mod¬ 
elled  errors  associated  with  dynamic  lags  in  the  gimbal 
angles  and  with  rate  gyro  transients  in  the  angle  rate 
measurements.  The  capability  to  turn  these  different  errors 
on  or  off  in  an  error  budget  gives  the  filter  designer  a 
powerful  tool  for  locating  the  areas  in  which  filter  design 


improvements  can  achieve  the  greatest  gains  in  filter 
performance  . 

The  filter  model  should  also  be  critically  reexamined 
during  advanced  testing.  The  angle  rate  measurement  equa¬ 
tions  in  the  filter  model  are  not  validated  by  performance 
results  in  this  study.  This  could  possibly  stem  from  inap¬ 
propriate  angle  rate  equations  in  the  filter  model,  mis¬ 
aligned  geometry  between  the  truth  model  and  the  filter 
model,  or  a  combination  of  both.  Another  recommended  im¬ 
provement  to  the  filter  is  to  model  the  dynamic  lags  in  the 
gimbal  angles  by  augmenting  extra  states  to  the  basic 
9-state  dynamics  model.  The  dynamics  equations  for  the 
extra  states  are  based  on  the  truth  model  equations  for 
calculating  the  lags. 

Finally,  the  advanced  testing  program  should  expand  the 
filter  tuning  procedures  followed  in  this  study.  For  exam¬ 
ple,  additional  noise  terms  could  be  added  to  those  dynamics 
equations  requiring  ownship  velocity  terms.  This  noise 
would  account  for  the  uncertainty  in  ownship  velocity  gene¬ 
rated  by  coordinate  t r ans f orman t ions  using  noise-corrupted 
INS  angles.  Another  possibility  is  ad  hoc  tuning  using  non¬ 
zero  off  diagonal  terms  in  the  driving  noise  or  measurement 
noise  matrices.  A  third  recommendation  in  this  area  is  to 
investigate  adaptive  tuning  techniques  again.  The  method 
used  in  Chapter  VI  to  increase  gimbal  angle  measurement 
noise  strength  during  ownship  maneuver  can  be  adapted  to  use 
a  variable-strength  noise  value  calculated  as  a  function  of 


range  dependent  noise  (glint)  and  roll  rate  dependent  noise 
(angle  lags). 

In  addition  to  continued  testing,  several  recommenda¬ 
tions  can  be  made  concerning  filter  implementation.  These 
recommendations,  based  on  the  evaluation  results  obtained  in 
this  study  are  as  follows: 

(1)  Minimize  complexity  by  implementing  a  9  to  12 
state  filter  in  the  fire  control  computer. 

(2)  Preserve  linearity  in  the  filter  dynamics  model  to 
avoid  the  computational  burden  of  on-line 
numerical  integration. 

(3)  Avoid  observability  problems  by  choosing  a 
coordinate  frame  in  which  all  modes  will  be 
excited  regardless  of  the  trajectory  flown.  The 
north,  east,  and  down  frame  used  in  this  study  may 
not  be  adequate  in  this  respect  under  some 
conditions. 

(4)  Implement  the  final  design  using  a  square  root  or 
U-D  algorithm.  The  impact  of  short  wordlength  is 
not  specifically  evaluated  in  this  study,  but  is  a 
well  known  cause  of  numerical  instabilities  in 
conventional  Kalman  algorithms. 

(5)  Tune  the  final  filter  for  several  different 
trajectories  to  insure  it  performs  adequately  for 
the  variable  conditions  of  the  air-to-air  combat 


environment . 


The  above  recommendations  are  not  all  inclusive.  Many 
techniques  are  available  to  today's  optimal  filter  designer 
and  all  cannot  be  examined  in  a  work  of  this  scope.  This 
study  has  demonstrated  the  feasibility  of  a  Kalman  filter 
for  the  F-AE/G  LRI  function  and  identified  several  critical 
issues  which  will  affect  the  success  of  any  Kalman  filter 
employed  for  that  purpose.  The  models  provided  in  this 
study  are  a  suitable  basis  for  experimental  filters  which 
can  be  further  improved  by  augmented  models  and  continued 
testing . 


Appendix  A :  Binary  Calculations 

Maximum  relative  range:  182,400  feet 

In  binary:  101100100010000000 
8  foot  precision:  1 0 1  I  00 1 000 1 000 - (14  digits  required) 

Maximum  total  range:  4,152,640  feet 

In  binary:  1 1 1 1 1 1 0  I  0  I  1  1 0 1 0 1 000000 
8  ft  precision:  111111010111010100 -  (18  digits  required) 

Maximum  relative  velocity:  1520  fps 

In  binary:  10111110000 

4  fps  precision:  10111110 -  (8  digits  required) 

Maximum  total  velocity:  760  fps 
In  binary;  1011111000 


4  fps  precision;  1011111 - (7  digits  required) 
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Q  (t  )  is  calculated  by  evaluating  (4-16),  with 
d  i 


Q(t  )  = 


‘N 
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0  0 


Carrying  out  the  indicated  matrix  multiplications  in 
(4-16)  and  integrating  over  the  0.04  second  update  interval 
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The  same  procedures  apply  for  the  GMA  model.  For 
example , 
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Evaluating  at  T  =  0.5  and  At  =  0.04  yields  the  values  listed 
in  (4-25). 
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These  coordinate  references  can  related  by  standard 


order  dependent,  Euler  angle  rotations.  For  instance,  t 


full  transformation  from  geographic  to  antenna  coordinates 
is  represented  by: 


=  [\  ]  [^-2]  [^]  [^]  [\] 

e  a 


where  [  A]  are  the  following  n-axis  rotation  matrices 
(38:20) : 
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where  superscript  T  indicates  the  transpose  of  the  above 
rotation  matrices.  Specific  transformations  used  in  this 
study  are  calculated  below. 


TRANSFORMATION  1:  F-4E/G  velocity  in  radar  reference 
coordinates  can  be  calculated  as 


(C-6) 


(C-7) 


(C-8) 


TRANSFORMATION  2:  Converting  a  vector  from  radar 
reference  coordinates  to  geographic  coordinates  requires 

Vp  =  [T21  Vp  (C-9) 

where  the  transformation  [T2]  is  given  by 


[T2]  = 


[^4-]'^  [^0]'^  [^-2]'^ 


(C-10) 


Carrying  out  the  indicated  matrix  multiplications  yields  the 
following  matrix  elements: 


=  {cos  'V  cos  6  cos  -2  }-  (cos  4*  sin  0  cos  ([)  sin  -2) 
-{sin4^sin<{)sin-2}  (C-11) 

=  {sin  4'  cos  6  cos  -2  }-  (sin  4*  sin  0  cos  (J)  sin  -2} 

+  (cos  4*  sin  (J)  sin  -2  }  (C-12) 

[2^^  =  {-sin  6  cos  -2}  -  {cos  9  cos  (1>  sin  -2)  (C-13) 

r2j2  =  {cos  4*  sin  0  sin  ({)}  -  {sin  4*  cos  (|>}  (C-14) 

r222  =  {sin  4”  sin  0  sin  <t) }  +  {cos  4*  cos  (J)}  (C-15) 

r222  =  cos  0  sin  ({)  (C-16) 


r2j^2  =  {cos  cos  0  sin  -2}  +  (cos  4^  sin  0  cos  (j)  cos  -2} 


+  {sin  4*  sin  <i>  cos  -2}  (C-17) 

^^23  ~  {sin  4*  cos  0  sin  -2}  +  {sin  4*  sin  0  cos  (t  cos  -2} 

-  (cos  4*  sin  (J)  cos  -2}  (C-18) 

T233  =  {-sin  0  sin  -2}  +  {cos  9  cos  (J)  cos  -2}  (C-19) 

One  possible  simplification  is  the  approximation 

cos  -2=1  (C-20) 


TRANSFORMATION  3:  The  inverse  of  transformation  2 
above  is  given  by 
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where  [T3]  =  [T2]  due  to  the  orthonormality  of  all  rotation 


matrices . 


TRANSFORMATION  4:  The  full  transformation  from  geo¬ 


graphic  coordinates  to  antenna  coordinates  expressed  in 
(C-1)  can  be  obtained  by  applying  two  additional  rotations 
to  the  result  of  transformation  3  above: 
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The  elements  of  this  matrix  are: 
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FIGURE  E.22:  GAUSS  MARKOV  MODEL,  STATE  1,  5  RUNS 

BEAM  TRAJECTORY,  0=5000,  ANGLE  RATE  MEAS.  INHIBITED  (ORIGINAL  R) 


163 


FIG 

BERi 


bill  L  UriTL  HNO  TIML  -  Bb/ll/lb.  lb. WPHtb,  DISSPLH,  SUl  LPl.  VLBblUN 


FIGURE  E.24:  GAUSS  MARKOV  MODEL,  STATE  3,  5  RUNS 

BEAM  TRAJECTORY,  0=5000,  ANGLE  RATE  MEAS.  INHIBITED  (ORIGINAL  RJ 
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FIGURE  F.l:  GAUSS  MARKOV  MODEL,  STATE  1,  5  RUNS 

BERM  TRAJECTORY,  0=1000,  ANGLE  RATE  MEflS.  INHIBITED  (ORIGINAL  R1 


FIGURE  P.2:  GAUSS  MARKOV  MODEL,  STATE  2,  5  RUNS 
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FIGURE  F.4:  GHU5S  MRRKOV  MODEL,  STATE  4,  5  RUNS 

BERM  TRAJECTORY,  Q-1000,  ANGLE  RATE  MEAS.  INHIBITED  (ORIGINAL  R) 


bU^  L 


E  r.5:  GAUSS  MARKOV  MODEL,  STATE  5,  5  RUNS 

TRAJECTORY,  0-1000,  ANGLE  RATE  MEAS.  INHIBITED  (ORIGINAL  RI 


BEAM  TRAJECTORY.  0-1000,  ANGLE  RATE  MEAS.  INHIBITED  (ORIGINAL  R) 


.lib. 


E  P.9;  GRLISS  MARKOV  MODEL,  STATE  9,  5  RUNS 

TRRJEUTORY,  Q-IOOG,  ANGLE  RATE  MEHS.  INHIBITED  I  ORIGINAL  R) 


AD-A164  2fS 
UNCLASSIFIED 


KRLHmi  FILTER  DESIGN  FOR  TIC  LONG  RANGE  INTERCEPT 
FUNCTION  OF  THE  F-4E/0  F. .  <U>  AIR  FORCE  INST  OF  TECH 
URIGHT-PATTERSON  AFB  OH  SCHOOL  OF  ENGI. .  R  C  HALBERT 
DEC  B3  RFIT/'0E/EN0/8SD-2S  F/B  17/9 


3/4 


MICROCOPY  RESOLUTION  TEST  CHART 
ciuR'Aj'  at  ‘•’•NOAorys  a 


FIGURE  P.IO:  GAUSS  MARKOV  MODEL,  STATE  1,  5  RUNS 

BEAM  TRAJECTORY,  0-100,  ANGLE  RATE  MEAS.  INHIBITED  [ORIGINAL  Rl 


. yuf'''  -f-.r-  f  i  'n  f.i  iwjh  -  *jii  ciuh 


FIGURE  F.12:  GAUSS  MARKOV  MODEL,  STATE  3,  5  RUNS 
BEAM  TRAJECTORY,  Q-lOO,  ANGLE  RATE  MEAS.  INAIBITEI 


FIGURE  F.13:  GAUSS  MARKOV  MODEL,  STATE  4,  5  RUNS 

BEAM  TRAJECTORY,  Q-lOO,  ANGLE  RATE  MEflS.  INHIBITED  (ORIGINAL  R) 
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19.  Abstract  (continued) 


Filter  tuning,  error  budgets,  adaptive  techniques,  and  observability  issues  are 
addressed  during  filter  evaluation.  Results  of  the  evaluation  indicate  the  filter 
designs  can  meet  the  requirements  of  the  F-4E/G  fire  control  system.  Recommendations 
are  made  for  continued  testing  and  for  operational  implementation. 
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